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TWO- BODY ERROR ANALYSIS PROGRAM 
Daniel P. Muhonen 


ABSTRACT 

The Two- Body Error Analysis Program is a rapid computer 
program designed to evaluate the resulting state vector and 
covariance matrix for an orbiter after one coast and one burn 
maneuver. After the burn, histograms are computed for the 
orbital elements and other state dependent variables, includ- 
ing delta-velocity corrections necessary to attain a prescribed 
orbit. 

The program is coded in FORTFAN IV for the IBM 360/95 
computer. This document contains both a user’s guide and a 
complete description of the analysis. 
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SUMMARY 


Initial conditions for this program include a state vector and 
covariance matrix for an orbiter, together with constant thrust 
vector and thrusting errors. The program will compute the 
resulting state vector and covariance matrix after one coast 
and one burn maneuver. In addition histograms are computed 
after the burn for orbital elements and othe/ state-dependent 
variables, including ae n a- velocity corrections necessary to 
attain a prescribed orbit. A two-body central force field is 
assumed throughout. To increase the flexibility of the program, 
an input option is included which will divide the burn time into 
a number of sub-intervals. This will reduce numerical trunca- 
tion errors for long burn maneuvers, and it enables the user 
to program variable thrust vector and thrusting errors as step 
functiors through the burn. The histograms are printed on the 
output printer in both graphical and numerical form. 
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TWO- BODY ERROR ANALYSIS PROGRAM 


INTRODUCTION 

The two-body error analysis program propagates a state vector and covariance 
matrix for an orbiter through coast and bum maneuvers. Histograms of the 
final orbital elements and other state-dependent variables are evaluated by a 
Monte-C'arlo procedure. Provision is included for calculating a histogram of 
the delta-velocity required to effect prescribed semi-major axis, eccentricity, 
and inclination (optional). 

RESTRICTIONS 


Installation 

System and Configuration 
Source Language 
Subroutines Required 


Functions Required 


GSFC 

IBM 360/95 
FORTRAN IV 

MTXPR, GAUSS, RANDU, EIGEN, 
BURNST, PARTAL, POWERX, 
MTRPLY, TWOBDY, CONVET, 

MTRX, TAB1, CROSS, ORB, TCONIC, 
HISTO, DELVS 

NS AMP, BARN1, DOT, ARKTNS, 
FNORM 


Storage Required 
Common Storage Required 
Input Range and Limitations 

Output Range, Limitations, and - 
Accuracy 

Running Time 


400 K BYTES TOTAL 
None 

Monte-Carlo sample size £ 10,000, 
number of histogram intervals < 50 

Calculations in double precision; two- 
body force field; linear error propa- 
gation; normal distribution of the 
Cartesian state vector; burn errors 
uncorrelated. 

0.25 minutes should be allowed per 
1000 Monte Carlo sample vectors 
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INPUTS/OUTPUTS 


Remarks 

All input is in card form described in Table 1. Initial state vector, covariance 
matrix, coast time (can be zero), bum time (can be zero), weight, weight loss, 
thrust vector, thrust errors, and histogram information can be input. A sample 
input case is given in Figure 1, and the corresponding printer output is included 
under Figure 2. All output is in units of kilometers, kilograms, seconds, and 
degrees, unless stated otherwise. Histograms of the following state-dependent 
parameters, after burn, can be printed: 

1. Eccentricity 

2. Semi-major axis 

3. Inclination 

4. Longitude of ascending node 

5. Argument of periapsis 

6. True anomaly 

7. Periapsis radius 

8. Apoapsis radius 

9. Delta velocity to attain a given semi-major axis, eccentricity, and 
inclination (optional). 


IMPLEMENTATION 

Remarks 

The program is available on tape with the source and the object on separate 
files. 

I/O Logical Unit Assignments 


FORTRAN 


Logical Unit 

Device Type 

Use 

Space Required 

5 

Reader 

Input 

— 

6 

Printer 

Output 



Control Cards or Language 

Included with the sample input (see Figure 1) 
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Table 1 

Input Card Formats 


Any number of input cases can be stacked; each case must include the following 
cards: 

Card No. 1 . . . Title for output, alphanumeric, columns 1-80. 

Card No. 2 . . . Namelist name, & T 13 DATA , columns 2-8. 


The remaining data is in namelist format, i.e., Parameter = ............. . 

columns 2-72. The following parameters can be input; if they arc not, either 
default values or previous stacked-case values are assumed. 


Parameter 

■ 

Dimension 

Default Value 

Description 

DTCI 

1 

0. 

Coas* time interval (sec.) 

DTBI 

1 

0. 

Total burn time interval (sec.) 

NI3URN 

1 

i 

Number of computation incre- 
ments in burn interval (see 
method). 

WI 

1 

0. 

Initial weight (lb.) 

WCI 

1 

0. 

Total weight increase during 
burn (lb.) (generally negative) 

PITCHI 

1 

0. 

Initial nominal pitch during burn 
(deg.) 

YAWI 

1 

0. 

Initial nominal yaw during burn 
(deg.) 

THRUSI 

1 

0. 

Nominal thrust during burn (lb.) 

ZDELT 

1 

0. 

Uncertainty for Monte- Carlo 
sample size (used only if 
1 < NCONF < 6) 

NCONF 

1 

0 

= 0 implies no histograms 
generated 

= 1 implies confidence level = .9 
= 2 implies t yiriidence level - 
.95 

= 3 implies confidence level = 
.98 
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Table 1 (continued) 


f 1 

Parameter 

1 

Dimension 

Default Value 

Description 

1 

NCONF 

(cont.) 

1 

0 

= 4 implies confidence level = 
.99 

= 5 implies confidence level = 
.995 

= G implies confidence level = 
.998 

> 6 implies sample size = 
NCONF 

AS 

1 

0 . 

= 0 implies no delta-velocity 
histogram computed 
> 0 implies delta- velocity 
histogram computed to 
attain semi-major axis = AS. 

ES 

1 

0 . 

Delta -velocity histogram com- 
puted to attain eccentricity = 
ES. 

| A IS 

i 

1 

-1. 

< 0 implies delta- velocity 

histogram includes no incli- 
nation changes, 

_> 0 implies delta-velocity 
computed to attain inclina- 
tion = AIS. 

SIGI3I 

3 

o 

«• 

o 

• 

p 

Thrust error a T (lb.), 
Pitch error a a (deg.), 
Yaw error a (deg.) 

IHIST 

6 

1, 2, 3, 4, 5, 6 

Flags to indicate which histo- 
grams are desired. 

Computations are provided for 
nine state dependent param- 
eters, and up to six of these 
can be chosen for histograms. 
The i™ value of IHIST im- 
plies that the i™ histogram 
will be as follows: IHIST (i) 
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Tabic* 1 (continued) 


Parameter 

I - j 

— 

Dimension 

Default Value 

__ _____ ^ „ 

Description 



I HIST 
(cont.) 

* 

i 

6 

1, 2, 3, 4, 5, (3 

= 0 . . . not printed 
= 1 . . . eccentricity 
= 2 . . . semi- major axis 
= 3 . . . inclination 
= 4 . . . longitude of ascending 
node 

= 5 . . . argument of periapsis 
= G . . . true anomaly 
= 7 . . . periapsis radius 
= 8 . . . apoapsis radius 
= 9 . . . delta velocity 

ZUB 

(3,6) 

ZUB (1 , i) =0. 
ZUB(2,i) =50. 
ZUB(3,i) = 0. 
i = 1, ...» 6 

ZUB (1, i) = beginning value of 
i™ histogram 

ZUB (2, i) = number of intervals 
in i™ histogram 

ZUB (3, i) = final value of i™ 
histogram 

(ZUB (1, i) = ZUB (3, i) = 0] 
implies that the limits for 
the i™ histogram are those 
of the sample generated. 

XII 

6 

o.,o.,o.,o.,o.,o. 

Initial state vector, j , in 

inertial cartesian coordinates 
(km., km./ sec.) 

PII 

21 

o 

• 

o 

• 

• 

• 

o 

Upper half of initial covariance 
matrix, input row -wise 
[e.g., PII (6) = Pj 6 , PII (7) = 
p 22 i. 

IPCOOR 

1 

1 

Coordinates and units flag for 
PH. 

= 1 implies local tangent coor- 
dinates, [(r x r)x r, r x r, 
r J , and units in ft. 2 and 
ft. 2 /sec. 2 
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Table 1 (continued) 


1 

Parameter 

Dimension 

Default Value 

Description 

IPCOOR 

(eont.) 

i 

1 

= 2 implies local tangent co- 
ordinates and units in km. 2 
and km. 2 /sec. 2 . 

- 3 implies inertial coordinates 
of XII and units in ft. 2 and 
ft. 2 / sec. 2 . 

= 4 implies inertial coordinates 
of XII and units in km. 2 and 
km. 2 /sec. 2 . 

XMU 

1 

398603.2 
( ij . of earth) 

Gravitional constant times mass 
of central body (km.Vsec. 2 ). 


FINAL CARD .... namelist data end flag, &LND, columns 2-5. 
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0.12749537074640-03 
0.22019699341500-04 
— O . 30999497099930-05 
0.30939497099990-05 


0.96494779965640-21 
0.19293499039490-21 
-0.71150767969960-19 
-0. 30939497099950-05 
0.34067109378660-03 
0.296461 531 53900-20 


-0.12)66366433390-16 
-0.29866397326710-19 
0.13501240169990-01 
0.30939497099950-05 
0. 1474369491 3450-20 
0.34062103578660-03 


COVARIANCE MATRIX (Lite TAM — F T 1 

-0.3853903935 101 11-13 -0.20999682441940-13 -0.79063409639290-10 0.91906422686440-14 

-0.20359074781190-12 -0. 11(166992274130-12 0.17763115441040-10 0.37796207293080-14 


0.60767439446290-14 -0.13741607742140-11 
0.1641673107492 0-14 -0.27664U126676D-12 


Figure 2. Sample Printer Output 



0.60*9*05121*3*0-10 0.675939069*0650-10 

-0.222773209(10970-13 -0.99800315969*10-1* 
0. *6990? 9*9722*0— l* 0.25303092**3*50-1* 
— 0. 111283979 76690— 1 1 -0. 13799109211020-13 

COVAKlANCfc MCTRIX I INERTIAL) 

0.1111 9628726890-02 -0.360399*722808 0-01 
-0.360399*7228080-01 0. 1 168229612*970 Ol 

0.21937697931690-08 -0.698229272*0610-07 
0.1 619*91 H337290-0* -0.92371028*67390-03 
-0. *1996*83782*20-03 0. 13*89089663990-01 

0.18098088298610-09 -0.989*21760*0770-0* 
EIGENVALUES 

0.1169*97098*9*0 01 0. 3*069 10789*9*0-03 

(1RHITAL ELEMENTS NOMINAL 

ECC SNA 

0.7316279*2912*0 00 0.2**968977167*0 09 

TRANSIT1L MATRIX 

0.63939699237190 01 -0. *3908609019790 02 
0.192371*1986190 01 -0.96960681313**0 02 
0.70119*92671360 01 -0.2**793*6970000 O? 
0.31193*362*00*0-0* -0. 13380?97700?2'’-02 
0. *709*763332 390-03 -0.72919806779010-0? 
0.23099*89988610-0* -0.8*929013281730-03 


AFTER COAST 


STATF VECTOR 

0.1*3731962*7700 0* -0.*2329069899760 09 

COVARIANCE MATRIX (LOCAL TANGENT) 
0.12**7010162380 06 -0.26010107739690 01 
-0.26010107739900 01 0.727696*6229860-01 

-0.916790022*9030 09 0.2808*911100200 01 

0.31*61691861260 01 -0.97*901*1278710-0* 
-0.27979559979980-01 0.77196319959270-03 

-0.13270031138*90 02 0.363*67693016*0-03 

COVARIANCE MATRIX HOC TAN — FT) 
0.1339785*855930 13 -0.27997056621380 Oh 
-0.27997056621650 08 0.7832*316796730 06 

-0.98682*90766210 12 0.302303971*2290 08 

0.3386509*956970 08 -0.10*89**9526510 O* 
-0.29686396799820 06 0.H3050399829970 0* 

-0.1*2837*7567310 09 0.39123350369*90 O* 

COVARIANCE MATRIX I INERTIAL 1 
0.9060*0*6*9*750 05 0.80676763*71960 05 

0.80676763*71960 05 0.9*920195215890 05 

0.50577512216680 05 0.65*21677139990 05 

0.1963119899*300 01 0. 23920 1 1 H 697760 01 

0.1132*06090*110 02 0.12368239662130 02 

0.12295783752930 01 0.1*72090*9*5*70 Ol 

EIGENVALUES 

0.20023300*9*020 06 0.13531379396690 05 

ORBITAL ELEMENTS NOMINAL 

ECC SMA 

0.7316279*2512*0 00 0.2**568977167*4' 05 


AFTER BURN 


STATE VECTOR 

0.15399929356890 0* -0. *2320203278200 05 

COVARIANCE MATRIX (LUCAL TANGENT) 
0.98*905313186*0 05 0.512985905*6810 05 

0.512985905*6810 05 0.26T1899235181O 05 

-0.819*2655176680 05 -O. *266059*136*90 05 
0.2*998056*71500 01 0.130080838*6820 01 

0.12659319059560 01 0.6602158*212030 Oo 

-0.118*706*369720 02 -0.616868*2901**0 01 

COVARIANCE MATRIX (LUC TAN — FT) 
0.10601*36217020 13 0.55217362362050 12 

0.55217362362050 12 0. 2876009393071)0 12 

-0.88202370388150 12 -0.65919*972007 10 12 
0.26907693326870 08 0.1*00178975*060 08 

0.13626382329750 08 0.71065066316510 07 

-0.1275207835667D 09 -0.66399187961790 08 

COVARIANCE MATRIX (INERTIAL) 
0.907913008*9070 05 0.81335669316930 05 

0.81335669316930 05 0.961125*0162*20 OS 

0.50608973226790 05 0.657979616*9580 05 

o.l 9*| 20** 7751 20 01 0.2383957*250510 01 

0.113JTT070601890 02 0.12*93036**1780 02 

0.121|83193136730 01 0.1*696053*63920 01 

EIGENVALUES 

0.2016*991973500 06 0.13560210250670 OS 


BURN ERRORS 


STATE VECTOR 

0.15399929356890 0* -0.62320203278200 05 

COVARIANCE MATRIX (LOCAL TANGENT) 
0.1006910*50*380-01 0.118985513*2270-02 


n. i?s«h6R«ooooon or o. I3?nn3660unnou o* 
O. ) 3 ?(j03 66000000 O* O. 23697725000000 03 
-0.0**9»3(l7>6*n?P-l? —0.33303009000000 02 
0.1*532619(100000 06 0.33303009000001) 02 


0.2153765793 1690-0R 0.1615*51*337260-0* 

-O. 69*^2*272*06 in -07 -0.5237102**67390-03 
O.* I 73 139 6**9080—1* 0.31 301017662630-10 

0.3 1301017662530-10 0.92*9*060**1660-0* 

-0.80597*05 767960-09 -o. 1037209*067 160-0* 
0.3*9 M 9377*10110-11 -0.1 31*96267*7260-03 

O. l**72R067212*n-03 O.219567i>095«7?D-o* 


INC UMFG* 

0.7*5280060017*0 02 0.27176679*5*7*0 03 


0.6?5051H»*5?75n 01 0. *275*3*5*003*0 06 

0.205*1076*1**50 00 0.56063*6*157120 06 

-0.3052650956*830 01 O . 2* l 772 /9 | 16 ion 05 
-0.72)166*71636*0-05 0.13*6328910*170 01 

0. M3fc0**39?*?on-03 0.7121*067167330 ol 

0. 151u911*36?*8n-05 0.95679166*00650 00 


0.71 585*5*568*90 02 0.1395007203**60 01 


-0.9)6790022*5030 05 (.. 31**169 1 *6 ) 2*0 01 

0.2*08*911100*30 01 -0.97*501*17 7*960-0* 
O. *929*7*0*97320 06 -0.30*75710*96*30 01 
-0.30*75710896*30 (>1 0.1067*59*727*20-03 

O. 79779552*10850-01 -O. 103330? 750* 200-05 
o. 11877*770792*0 0? -0.*U9*55516 u7*5D-03 


-0.9*682*90766710 1? 0.33*6*09*966970 08 

0. 302303571*25 30 0* -0.10**9**9*265*0 ()* 
0.9611559675*650 17 -0.3323*35007*690 08 
—0.3323*3 5007*690 O* 0.11*9*3*7*671*0 0* 
0.3206**5**97130 06 -0.11127382079*10 02 
o. 127H67***°?*50 09 -0. **116*95*06100 0* 


0.505775122166*0 05 0.1963119*99*3181 Ol 
0. *5*21677139950 05 0.2392011*697760 01 
0.2*2*0213176090 05 0.11062093507320 01 
0.11062093507320 01 0.62370*12*00020-0* 
0.63595635992910 01 0.3091*977*70*60-03 
0.693*5*2670*130 00 0.337*0313709*60-0* 

0.72682*5*3*3360-01 0.12*665021613*0-07 


INC UMFGA 

0.2*52*0060017*0 02 0.271 766795**7*0 03 


0.929*326*967890 0? 0.3074136***7090 01 


-0.* 19*2653 l 766*0 05 0.2*99*056*71500 PI 

-O. *>66039*136*90 05 0.130O*0M3R**8?O 01 

0.9(100076996*030 05 -0.277075 733616*0 01 
—O. 277075 73*61 6*0 01 0.870763623296011-0* 

-O. 1*0175631 10740 M 0.39816853678030-0* 
0.11961885060670 0? -0.36732327807250-03 


-0.88202370388150 1? 0.26907693326870 OH 

-0.63919*97200710 12 0.1*00178975*060 08 

0.968760559*9990 1? -0.29*2*19*0*0760 OH 
-0.2982*19*0*0750 08 0.9372*237920330 03 

-0.130*838*363170 08 0. *2857**2099100 03 

0.12*756703622*0 09 -0.39538362230120 0* 


0.506*8973226790 03 0.19*120**773120 01 
O. *57479616*9580 05 0.2383957*250510 01 
0.28306*52*23010 04 0.109391257*0310 Ol 
0.109391257*0310 01 0.6126576*977200-0* 
0.6390235*98*680 01 0.30712227*29*30-03 
0.68766563065630 00 0.3311*200232320-0* 

0.16539223780770 00 0.1236539225*310-07 


0.929*3268967690 02 0.307*136*6*7090 Ol 


-0.252367213819*0-0* 0. *7305*06 179930-03 


-l>.765R*l)75?37«9|)-12 0.1*532619000001) 06 

-0.3 330 30090(10000 02 O . 33303(l»l9O(l!)PUO 0? 
0.3666 4 155*00001) (I* O. 1 58699*7x98 79’)-l 3 
0.31910*6*682*50-13 0.3666*165*00001) 0* 


-0.*1 506*8 3782*2(1-03 0. 1 HO 58(1* *2 986 10-0 6 

0. 1 3**60*6463560-01 -0.585*21760*0770-0* 
-d. *0697*06767960-09 0.3*9*9377*10110-1 1 

-0.1037209*067160-0* -0.131*96267*72>o-03 
O. 3*0? 172*636321)— 03 -0.55*6*6) 3*843*0-05 
-0.56*6*51 3*853*0-05 0.2705*666511*311-03 

0.6)936107099**0-15 U.U 

AK GP TMFTA 

O. 1*1100000051120 03 0.65370*56)6601111-05 


-if. 7061*) 36651660 O* 0.2311“3?7317*ll) 06 
0.1*393223622760 (1* 0.3U*863*O6*12,»l) IIS 

-H. 3* ) 231 26960750 ()* 0.1 30577*0191* (.1 11* 

(1.5721 307**5*770-01 11.83*39*1 661 9660 0(1 

-o. 1 1*178*1 5**951, (10 (1.3*6996067*3311) HI 

0.31749*0369*160-01 0.368*07**2**7*11 (10 


0.6*000813657*70-01 0 . 7890?? 36767720 OP 


-0.2767966567998(1-01 -II. 1 3270031 1 3**91. 
0.77166 319 459260-1)3 O. 363*6 76U3p 1 jpi-iO 
(1.29779 44?* 1 l»8 40— (» 1 0.11*77*770792811 02 

-o. 10 3 33027 4(J 8 2(1 0—05 -0. *09*656) *o /*•,' 1-1)3 
O. *1*1 19267*7*90-05 0. 3*53*19 30* 8) 33ii-ns 

(1.3*4 399 30*8 1 33|/- 05 0.161*0761 12116*0-01 


-0.296*639*799*211 (l* -O. 1 *2 * 37* 76*7 31 •) 09 
0. *31150399*299*0 0* 0.391233403*96/11 II* 

0.3204**5**971 3|l 118 o. 127*s?**huv*iii pu 
-O. 111223*2079810 02 -O. **1164 <16 8.18100 n* 
0. **06144*3751(1, 02 I).* 1*8*04027*8 / i| i>7 

(I. *1 ***1151)726670 02 11.17*16*32 33 4/. /•! ..s 


o. 1 1 32*06091)* 11 O 02 il. |22947*37S?9.|| ol 
II. 123**230*821 3|» (I? 0.1*/?l’«O*9*s*/ll ,l| 

0.63494*34992910 HI P, 493*6*2670*1 SO IIP 
11.309) *977* 7u**U-iM <1. 337*0 3 1 37098*. i-o<. 
0.16*141 3568(1 1 (>||-P2 (I. 19(,U7** 734 ms /•(— (.8 

O. 19092**739887(8103 II .29 1 4 8 884.. 0*2. i *1—0* 

II. *917*1 2**113**0-10 0.0 


AMGP Tms|h 

0. 1*0000000 1*9 7 1, in -0.1797972)0733*1 .13 


0.1*92290*3*81*0 OU I). 527106*7782*3.1-01 


0.126493191)898*11 111 -0.1l8*7P6*3*u/2 l i>7 
0.6402) 4**?17H»0 OO —0.6 1 8 86**2 40 1*/.* i ..1 
-O. 1*0176631 107*0 Ol 0.1196)8*40*06/0 07 
0.39* | 4*53678030-04 -0.3673? 377M077SO-O3 
O. 2*24992 *7 30 7*0-0* -O. l*59l*7«*n6334i)-o8 
-0.18 491 *790* 3340-03 0. 162*1230790*811—02 


0.13*2*3*2 329 76.1 Oh -(). 1275707A344* 7U 09 
0.71064(1**31*410 oV — 0.66 39 u l *78*1 7 VO 08 
-0.140**3**4*3171) 08 o. 12874*703*22*11 09 
0. *2*47**2owloo 03 -0.395 3*3*2? 301 2 *) 0* 
0.30*1*7*?***) All 03 —I). 2(H) 11708*1*2 O') II* 
-0.2001170**1*200 0* 0.1752*977019*60 05 


0.11377070*01*90 02 0.1218319413*73.) ol 

0.12*9303***17*0 02 0.1*69*04346.3720 111 

O. *3Q0?34OO****o 01 0.6*76*4*30*4*3)) OO 

0.30712227*29*30-03 0.3311 *200232 32 0-0* 

0. 1*433***4174*0— ()2 U. l9(MI()7*?3«*6«,. t -03 
0.190007*23**6*0-03 0.2**069*1 3771* : 

0.89752801*179*0-10 0.0 


0.1*922*0*3**1*0 OO O. 52 7 1(1*477 *2*3')“01 


0. 44704*936*9030-0* -O. 1*507*3*** 3080-05 
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SVS/.OHF-hV 

KFPT 

|FF/»Bt 

VOL SF* mis* 0 1 SY Si*. 


IFF/*b I 

SYS/. I"'— Y 

KFPT 

IFF/MBl 

voL sfp «uis* r.isvy.. 
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KFPT 
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vol sfp nos* r.isvy*. 


IfF^hbI 
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IFF/MB 1 

SYS1 .TFttFL 1 * 

KFPT 

IFF2*bt 

VOL St* MIS* r.isvsr.. 
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METHOD OF ANALYSIS 


The following parameters used in the analysis are input variables. All other 
parameters are built-in constants, or are calculated from the input data: 

At c Coast interval 

At b Burn interval 

W ( , Initial weight 

AW Weight change during burn 

a Nominal pitch during burn 

Nominal yaw during burn 

T Nominal thrust during burn 

—4 

X Nominal initial state vector (X 0j ) 

P 0 Initial state covariance matrix 

o\_, ct u , a y Errors in thrust, pitch, yaw 

The initial state covariance matrix can be input in local tangent coordinates, but 
all computation is made in an inertial central body Cartesian framework. This 
coordinate system applies to the analysis presented. (See Figure 3.) 

Propagation of the State Vector and Covariance Matrix Through Coast 


This is accomplished by integrating the equations of motion for the two-body 
problem through the coast time interval. Let t 0 be the initial time. We define 
t , , Xj, P 1 to be the time, state vector, and state covariance matrix at the end 
of coast (t j = 1 0 + At c ) . 

The equations of motion are: 


• • 

r " - M r / r 3 , 

where r (t) represents the position vector at time t, and ^ is the gravitational 
constant times mass of central body. 
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A 



k = [ ' 



— • _ — ♦ < - 4 , 1 

• " f ( r x v) x r ] 

Figure 3. Local Tangent Coordinate Frame for the Input Covariance Matrix. 


(x, y, z) Represents the Inertial Coordinate Frame of the 
State Vector 

(i, j, k) Represents the Local Tangent Frame. 
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The integration is accomplished by making use of the subroutine TWOBDY pre- 
pared by W. H. Goodyear.* Goodyear defines a new variable ; by the differen- 
tial equation: 


i = i/r, 


where <p is zero when t is t 0 . Then r and r are expanded in a Taylor series 
about r (t 0 ) in terms of <P. From these series expansions it is possible to ob- 


d X 


evaluated, and the state transition matrix <1 >j = 


/? 

(t,)\ 


1 ( - 


* (X ,.) is 


<‘l)/ 


1»\ 



— 1 is 

found 

by computing 

0 } ! 




the partial derivatives of the terms in the Taylor series. 
The covariance matrix Pj is then evaluated by: 


P 


1 


I), P 0 * 


T 

1 


Propagation of the State Vector and Covariance Matrix Through Burn, wi thou t 
Addition of Burn Errors 

Let t 2 , X 2 , P 2 represent the time, state vector, and state covariance matrix 
(without burn errors) after burn (t 2 = tj + At b ). 

The equations of motion are: 

U _ - / 3 ^ 

r - “ m r r J +— , 

m 


where m(t) is the mass of the vehicle and X (t) represents the thrust vector, 
evaluated by: 


* Goodyear, Vt. H., "A General Method for the Computation of Cartesian Coordinates and Partial 
Derivatives of the Two-Body Problem/’ NASA CR-522. 
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—4 

X 



“4 V *4 

( r * r ) x r 

• • 

(r x r) x ? 


? x r 


r x r 


T cos \fj cos a\^ 
T cos \p sin a 
T sin 0 


/ 


An approximate solution to the equations uf motion is obtained by integrating: 

-t -* 

” M r j A \ 

r 3 m ' ’ 

1 l 


where Xj = X (t t ) , and r x = r (t x ). 

The solution for this equation through burn is: 



where r p m x are the position vector and mass before burn, Am is the mass 
change during burn, and r 2 is the position vector after burn. This provides 

-4 ‘ 

r 


the state vector X 2 = 


= (x 2l ). The state transition matrix ^ = 


3x 2i 


is found by computing the partial derivatives from the above equations, 


3x u. 

and the covariance matrix P 2 is evaluated as before: 


P 


2 


<D 2 Pj 4> 


T 

2 ’ 


For long burn maneuvers provision is made to reduce the error resulting from 
the approximate solution to the equations of motion. An input parameter is 
available which breaks up the burn interval into n b sub -intervals; the state 
vector and covariance matrix are then propagated as above through each burn 
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THRUST 


ORBIT 


; > o) 


’itch and Yaw Angles. 

'he Pitch, a, is Measured Toward-j in the (i, j) Plane (Orbit 
Mane) from i to the Projection, (\., v, 0) of the Thrust 
'ector, Yaw, \p, is Measured Toward k from (K, 0) 


19 


step. Variations in the thrust vector through burn can easily be programmed 
into this scheme, if step functions are assumed. The computation of P 2 pre- 
cedes as follows: 



"b 


IT 



p 


2 


<!> P 
* 2 1 1 



The Inc lusion of the 1 Burn Krrors into the Covariance Matrix P 2 

Using the input errors - T , cr a , a diagonal covariance matrix for burn errors 
is assumed: 


b = 


0 

0 


0 


0 



This is transformed into a thrust covariance matrix first in a local tangent 
Cartesian coordinate frame and then in the reference inertial coordinate system, 

3 x, 


2 i 


Finally the state-thrust transition matrix " J is used to evaluate the state 

covariance matrix due to burn, and this is added to P 2 to obtain the final state 
covariance matrix, P 7t after burn. Let 


cos a cos \p -T sin a cos i// -T cos a sin \p ^ 

U - ( -sin a cos \p -T cos a cos \p T sin a sin \p 

sin >// 0 T cos v// J 


T = 


Cr i x rj) x t 
• • 

(T, » r j ) « r, 


r i K r i 


r l x r i 
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and 


Q 


d x 


2 i 


a \ 


ij 


where x = (A.) is the thrust vector before burn. 
The equations of Q are: 


a r 

T\ 


i. \ \ i A t , 2 

2 1 1 i 1 j b 

l j \ 2 A m 


9 - 


2m j + Am 
m 


log 


m j + Am 


m 


♦ b 


A_t l 

‘ J Am 


m . + Am 
A m 


log 


mj + Am 


m 


- 1 


a i 

T\ 


2 i 


lj 


' li ' lj A **> 




Am 


A m 


m . + Am 

- log I 

m j + A m \ m j 


+ 8 


* j Am 


log 


m , + Am 


m 


where b i . •- 0 , i ^ j and 5 x . = 1 . 

The final covariance matrix after burn is computed as follows: 


P 2 = P 2 + Q T U B U T T T Q t • 


If the burn interval is divided into n,, sub-intervals by the input option, the 
equations for the covariance propagation become: 


e = 


"r, 1 m 


T TT <t 
icn » 

Q, T k u. 


+ Q T U 

v "b n b "b 


P 2 = p 2 + e b e 1 
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Computation of Histograms of State Dependent- Variables 


Let (X) represent the vector of state dependent-variables for which histo- 

— <♦ 

grams are desired. The distribution of (X 2 ) is computed as follows: 


a. The Eigenvalues (E,) and associated matrix of Eigenvectors, R, is first 
computed for I* 2 : 


ding (EJ = R-'PjR. 

b. A set of n normally distributed random vectors, {A/3. : i = 1, . . . , n • , 
is generated with mean (0) and variance (E t ). 

c. A reverse transformation is made on each random vector, and the 
mean is shifted to the nominal by addition of X 2 : 


y i - X 2 + R A /3 i , i - 1 n. 

d. Finally { f (> t ) : i = 1, . . . , n } is computed, and the distributions of the 
state dependent-variables are presented as histograms of the elements 
of the vectors {£ (y i )} . 

Computa tion of the Delta-Velocity required to effect a given Semi-Major Axis, 
Eccentricity and Inclination 

To change semi-major axis and eccentricity, thrusting is assumed to be at 
apoapsis and periapsis. To change inclination thrusting is assumed to be at one 
of the nodes. If no inclination is desired, h v is calculated for both apoapsis- 
periapsis and periapsis-apoapsis burn sequences. If an inclination change is 
desired, twelve possible burn sequences of apoapsis-periapsis-node are con- 
sidered. The minimum of all possible Sv calculations is accepted. The follow- 
ing identities are used to calculate the Sv’s in the burn sequences: 

r a = a ( 1 + e > 


r p = a ( 1 - e) 
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A 


r z 

r r 

p « 

1 ±| 

( -- A 

cos w 

i 

n 

a 

\ a ) 





— 


v, - v( i ) 
where: 

r t = Apoapsis radius 
r = Periapsis radius 
a = Semi- major axis 
c = Eccentricity 
v = Velocity 


u = Gravitational constant times mass of central body 
r = Radius 
r n = Radius at node 
w = Argument of periapsis 
v i = Velocity to change inclination 
b i = Desired inclination change 
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FLOWCHART 
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The following is a FORTRAN source listing of the entire program. 


The listing is arranged as follows: 

Page 

1. Main program 27 

2. Subroutine MTXPR 30 

3. Function NSAMP 30 

4. Function BARN1 31 

5. Subroutine GAUSS 31 

6 . Subroutine RANDU 31 

7. Subroutine EIGEN 31 

8 . Subroutine BURNST 34 

9. Subroutine PARTAL 34 

10. Subroutine POWERX 35 

11. Subroutine MTR PLY. 36 

12. Subroutine TWOBDY. 36 

13. Subroutine CONVET 39 

14. Subroutine MTRX 39 

15. Subroutine TAB1 33 

16. Function DOT 40 

17. Function ARKTNS 41 

18. Function FNORM 41 


LISTING (Continued) 


Page 


19. Subroutine CROSS -41 

20. Subroutine OR II 41 

21. Subroutine TCONIC 42 

22. Block data IDCAM 42 

23. Subroutine HI STO 42 

24. Subroutine DELVS 43 


C '/«!*» PtOW.R AM 

C 

L ThFRP -- T 1 - nOfiY mu* a ri vs|v w»nr,« 

C 

ILLICIT «f- Al** I a-h ,u-v » 

MFAl*A WAHpS 

i • si i < ■ . 1 . i o i # i i* i t i • i , • i ■ i 

I • | . I. ' . I • . • 

UI'P*SIH»' /PCT I * II , / ST/ T S ( A I , T * T1 I IO|,HT{Q| 

DIM 'S I on W(6,A),o<j, 31**16), Fw(A t A|,7H ln<)f <• ,7 | 

D|’ H Sin" 70*1 3, 6 ) , ZFk* •«•»<») 

OI"t- >s Ilf IV(3|,ZX(5n) , KM a* , AIM » 
l) | <• ts | f l»« /HUM 
DlwFUSIO* AhFSI6| 

D|M* 'SIOM VI6.6) 

0|Mp SIMM ((6,6) 

im fns i imi st i m , x 1 1 < a » , p i m ?i i »sir,«m I 

OAT A FIT /3HFCC* 3H5 1 < . $H| Mt , amui r r. M , <,HAPr.P ,SmT PF T A , 3MFF 9 , AH A<*ur. , 

2 4MUFLV/ 

DATA m AW*S/iKFCC, 3MS ■ A , '*M| N( ,swili ■ (• r,A,<,MAMf.w ,ShThM«/ 

.. II i I J i /THfUTA/OTCI fftT'lfbl fWCt «»ITr.HtfVAM|fTM p . > 01 i i , « <>• , , 
AS, Cl, MS, BIIN fSlCnj tlMBf » 1 1 1 l‘IST,P|! fXMHf |*CHn« 

Call f rr s*- t i ? in, 2*6,-i , i • i i 
DTC I • 0.00 
OTA | ■ 0.00 

wi • o.no 

mci ■ o.no 

PITCHl » 0,00 
YAW | « 0.00 

T HBOS I ■ 0.00 
ZflFLT - O. 

MCOAF « O 
as ■ o.no 
FS ■ 0,00 
I Pf.fK IB « l 
A IS ■ -1.00 
MHORN = I 
S I Oh 111) ■ 0.00 
SIG«I<2) - o.no 
s |G*| I 3 » » o.no 
00 *95 I - 1 . A 

zo*( l.n ■ o. 

Z0KI2, I ) "50, 

ZU8C 3, I I « 0. 

X III I ) « 0.00 
1H|ST ( 1 ) « I 

895 COOT IN* IF 

00 *99 I * 1,21 

pi I m * o.no 

899 CnwTlMUF 

FACTO* * .00AAA*22l7'>u/.o».9B*'AAAnu 
X-O « 39«AL3.20(< 

1U RF AD ( 6 , 5(*0, t- .1-99991 I T I T L F (I) , I « 1 * lO ) 

B F AO (5,T DATA) 

DTC - OTC I 
OT* « OTP I 

W s W| 

WC « «CI 
PITCH » PITCHl 
YAW » Y A u 1 
THRUST * THRUSI 
SlGttl 1 ) = SIGH! ( 1 ) 

S 1 Gtt ( 2 ) * S I GH I I 2 ) 

S I Gh ( 3 ) » SIGH (3) 

00 *97 1-1,6 

mm - xiiii i 

00 *96 J - 1,6 

L = m 13-1 I / 2* J-6 
PH I, J) - PII(L) 

P 1 1 J , I ) > P 1 11 , J ) 

896 CO- TINIIF 

897 CONTINUE 
NP * wcnwF 

IF (MP.LF.O) GO TU 15 
IF | Mp . L F .6 ) IP - i’SA-p I ^ COMF . 70clt j 
IF INF. GT. 10000) .\P = lOOOC 
15 CONTIAMIF 

WR I T F ( 6 , 600 ) (TITLF(I), 1-1,10) 

wR I T* (6,610) l>TC,l)Tl ,' ,wc,P| TCH, YAm,TwriiST ,ZOFLT ,.'Cl" c , I Hi ST 


loon 

?onu 

3000 
4000 
6 < MW i 
AOOO 
PINK* 

HOOO 
9 00(1 

) oooo 

I l MOO 
1 2000 
| 3000 
1 AOOO 
16 000 
1 6 1 no 

1 6 0OM 

17000 
1 “ 000 
19000 
20000 
21000 • 

21 100 

2 2000 
2 3000 
24000 
26 000 
26000 
270oo 
2*000 
29000 
30000 
3 | 000 
3110(1 

31200 **** 

32000 

31000 

34000 

36000 

36000 

37000 

3*000 

30 00.) 

4 0000 
41000 
4 2000 
A30OO 
440O(> 

46 00O 
4 AOOO 
4 7 000 

4 AOOO 
49 000 
60000 
6 J 00( 

6 2000 
6 3 000 
54000 
56 000 

5 6000 
57000 

5 *000 

69000 
60000 
61 ono 
62000 

6 3000 
64000 
65 non 
66000 
6 7000 

6 AOOO 

7 2000 
7 300(1 
74000 
75000 
76000 
77000 
7 A 000 
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t 


r> o non 


WAllfc (6,605) AS,fc S, Al S»* MlH , XMI 
Hi|Tf is^IM H/UMI.Jl* J«).M,I»1,3) 

WRlTfc 16,620) I SlUBl I ) . I ■ 1. 3 ) 

WR | TP 16,630) Mill | f I* If hi 

WRITE 16,640) IPCUUH.I <P|< I • J). J>1*6) »|>1.6) 

fMAURN ■ NMIIRN 

OTH ■ OTR/ENHURN 

WC ■ Wf /ENAUKN 

AIS • Al S/57. 2957745100 

CALL MPHlX I ,XI (4) ,X II, f ?) 

(Ml 1R I « 1*6 
XC I M ■ <1111 
XHI I I ■ X | (M 
no 17 J • 1*6 

IP i |PCnnn,LP.?l PiilfJi • >' I ( I * .))•. onov <.«•*/ 

PI 1* J) ■ pI 1 1 • J) 

VII , J) ■ O.IX) 
hi | » j) ■ 0.00 
17 CONTINUE 

hi i • n ■ 1*00 

1M COOTIMIP 

I) I IPCOOR.fO. UOR.IPCOnR.fC.il CAtt Cl Vf TIP I »ll »X I (4) • . • I i 

WHITP 16,65(0 

CALL MT XPK I X I ,P 1 ,fc ,PV ) 

HR | TP 16,7001 I H AMP SII)*I*1«6|,(C 21 I I* I ■ 1 , 6 I 

psj • o.no 

IP IllTC.EO. 0.001 (.0 TO 66 1 

CALL TWOKOYI X I ,OTL ,XHU ,PS| , XC , P, C 1 . C 2. f 3 *CA , f.6 ,C6 ,C7 I 
WHITE 16,6651 I I P I I , J), J- 1,6) ,1 . 1,61 
CALL MTRX |P,P I ,P,6,6, - I) 

WRITP 16,660) 

CALL MTXPMIXC.P.P.tV) 

CALL ORRIXC.XCI* ) ,XMO,C?> 

WPITP 16,700) I NAMES! I ) * I • 1 , 6) , I f t I I ) , | * 1 ,6 > 

§91 |l ( OTP. 1 1 . M ) 6U i 

Oil AA6 J« 1 , 6 
C2I J)«XCI j? 

AB6 CUNT INOP 

Sir.ttll) ■ | Sir.RI 1 )».(lOAA4H??l 700)**? 
s|gai?I « !SIG*(2l*.0174S129100t**2 
SIGHI3) « I 5< I G8 1 3 ) * .0 1 746 329 3 OO 1**2 


74000 • •** 

ROOOO 
A 1000 
R ? OOO 

R 3000 • *•• 

6 3100 

A ?200 
A 32 1 0 
A 4000 
R6 000 

• ***I)PL FT I ONS 

R 7 000 
A 8000 
A 9 000 
40000 

9 1 000 ***• 

9 2000 
4 300O 
9 3100 
94 000 

94100 

96 000 

96000 • **• 

9 7 000 
9 8000 

04 000 
100000 
101000 
111 2000 
103000 
1 1>4000 
1 06000 
106000 
I 07000 
108000 
1 09000 
110000 
1 1 looo 
1 1 2000 
1 1 3000 
1 1 4000 
1 1 6 000 


on 30 I ■ 1,3 
1)0 20 J ■ 1,3 
HI I, J) ■ 0.00 
20 CONTINUE 

HI I, 1 ) ■ SIGH I I ) 

30 CUNTINIIP 

0(1 MAR I " 1 «NBLIRN 

CALL HMRNSTIC2.C2IA). XH, XHI4), 4,0, 00, Wf, PITCH, TV, OTH, THRUST , YA'-' , 


2 X MU ) 

X*>ASS ■ W*P ACTOR 
XHOOT « WC/0TH*PACTUR 

CALL PnWFRX I C 2 * C 2 I 4 ) , TV, XMASS, XMOOT , 0.00, XMll ,1 > T H , O ,H ) 
CALL MTRX|0,H,H,6,6,01 
CALL M TRx (0,V,V,6,3,0) 

CALL PAPTAL I P AH , C2 » C2 I 4 ) , T 60, P | Tlh , Y A w , THHUS T ) 

CALL MTHPLY I T60,PAR ,U, 3, 3, 3, 6, 3,6) 

CALL MT P PL Y |R, II, T SO *6* 3, 3, 6, 6*6) 

DO 676 J « 1,6 
Of) 677 L « 1*3 
V I J * L ) ■ VI J,L)*T5U( J*L> 

677 CUNT I NUF 
676 CONTINUE 


INSFHT HIJHN COAHPCTIONS HFHP MhFM THPY'kF P P A OY . .. 

IP I YAWI .fcO.0.00) GU TO 57R 
VinV2»C2l4)*XR|4)4C2l5 )*XBI6)+C?I6)*XAI61 
AmSV 1 *OSORT I C2 1 4 ) ** 2>C? I 6 ) *<■ 2 + C.2 I 6 1 * * 2 ) 
AHSV2»nS0RT(XHl4)**24XBIS)**?4XR(6)**?) 

HPTA«OARCOS( V10V2/I AHSV 1* AmSV2 ) > 

YAWb VAW4APT A/.O 1745324 300 
57A CONTI NUF 


DO HA7 J= 1 , 6 
C2 I J ) =x A ( j » 


1 1 6000 
1 17000 

I ■ 

1 1 9 OOO 
120000 
121000 
12 2000 
1 23000 
I 41 
1 26 000 
126000 
1 27000 
12 AOOO 
1 7A10O 
179000 
1 30000 
1 31000 
1 32000 

132100 
1 32200 
1 3 2 300 
1 32400 
1 33000 
1 34000 
1 36000 
1 36000 
137000 
1 3*000 
1 30000 
140 (MV 
1 4 1 JOO 
1 42000 
1 4 

1 t 4000 
1 iSMOO 
’.46 000 
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887 Cl'NTIMIF 

W*W*WC 

MHH COMTINttF 

CAll MTRX |H,P,P,6,6,-1) 

CALL MT««|V,htV,M,-n 
JUlll' ( *>(6 70) 

CALC MI xPM (XP , P,P ,► V I 
WH | T P ( 6,6 mO) 

CALL MTXPK(Xb,V,F,>-V) 

no •>{» I • 1,6 

1)0 M) J ■ 1,6 
P ( I * J 1 • P ( I , J 1 ♦ V ( I , J I 
ad cunt |NIIF 
50 COM? I NOP 

9 R | T F (6,6901 

CALL P'TXPK(XM,P,P,PV1 

CAl L OHH|X>-, XP ( 6 1 , X • o, f?) 

PRlTP (6,7001 ( N AM t S ( l 1 . I ■ 1 » 6 I , ( f, * < | ), | ■ 1 , 6 1 
552 CDNT I MOP 

A I ■ C2 (3 1 /57. 2957796100 
APP • C2I5) /37. 295779* IOC 
C / ( h I - C/< 21* ( 1.0o*(/( 1 1 ) 

C 2 ( 7 1 • C2( 21* ( 1.00-0/ ( 1 1 I 

c? ( 9 1 ■ n,n<' 

IP ( AS, r.T, 0.001 CALL OPLVSf C 7| 2 1 , L2 ( A 1 , C ?( 7) ,AS,FS,A1 , A 1 S, APP , XM|j , 
19)1 

wRITF (6,7601 (PLTU 1, 1-7,91, (C?(l 1, 1*7,91 

IP (np.lp.o) r.n tu in 

00 200 J « 1,NP 
00 150 I ■ 1,6 
So ■ o.nn 

IP (Ell l.GT. 0.001 SU*USnBT| i ( | I 1 
All) ■ WARM! -1,-1, 12787, SO) 

150 COMT 1 NOP 

CALL Tk PL VIEV,A, XL, 6,6, 1,6,6,6) 

00 160 I ■ 1,6 

XC ( ! 1 ■ XM I I ♦XCI I 1 

160 CUM T I MOP 

(.ALL OMH(XC,XC(6),XH), A> 

A| - A (11/57.2967795100 
APP » 6(61/57.29577951 00 

A( 8 1 ■ A ( 2 1 * ( l.OO*A( 1 1 1 

A( 71 ■ A ( 2 1 * ( 1.00-aI 1 1 ) 

A I 9 1 ■ 0.00 

! MAS. OT .0.001 CALL OFL VS ( A ( ? ) , A ( 8 1 , A( 7 1 , AS , F S , A I , A I S, APP , XMO, A (9 I I 
00 180 I>1,6 
L » I M 1 ST ( I 1 

IP (L.LE.0.0R.L.GT.9) GO TO 180 
Z H ( J , I 1 - A(L> 

Z8 ( J, 7 1 ■ 1. 

180 COMTINOP 
200 CONTINUE 

00 250 I « 1,6 
L ■ I M I ST ( I 1 

IF ( L .LP.O.OK.L.GT .9 1 GO TO 250 
ZfcL T ( I 1 > C2( L 1 
N Z P - ZNP( 2, I 1 
NZX « NZP-i 

CALL TAAKZKI 1, I 1 ,Zb( 1 ,7 I, 1,70R(1, I ), 7FPP0,ZPCT,ZSTATS,NP, 1 ) 

IP ( ZUH ( 3 « I l.GT.Z Uh I 1,1 1 1 ftp TO 720 
ZIJ8* 1,1) « ZSTATSI 6) 

ZUH (3,1 1 » ZSTATSI5) 

220 ZDX - I ZOH( 3 , | 1 -ZU8( 1, I 1 >/(ZOH( 2,1 1-2. 1 
Z X ( 1 1 - ZUH ( 1, I 1 
on 230 J ■ 2, NZX 

ixui « zxij-ii+zox 

230 CIINT INOF 

CALL MI STO( TITLE, EL T( L) , Z PCT , ZX , Z STA TS I 2) , ZSTATSI 3) ,NZP,ZEIT (] ), 

1 NP) 

250 CONTINUE 

GO TO 10 
9999 CONTINUE 
STOP 

600 FOR'- AT ( 1 0 AH ) 

600 FORMAT ( 1H , ,6X , 10AH 1 

605 FOR ( /23X,2HAS,23X,2HFS,?2X,3HMlS,?OX,5HMbMRN,22X,3HXMU/ 

2 3025. 1 6 , 1 ?5,n7S. 16) 

61L FORMAT ( /22X,3H0TC* 22X, 3H0TH, 26X , 1 H« , 2 3 X , 2 hwC , 20X ,6 HP ITCH/ 

1 502 5. 16//22X , 1HYAW, 19 X ,6HTHRllST, 20x, 6HZ0PLT , 70X ,5HMC0NF , 


169000 

1 6001)0 
16 loon 

151100 
161200 
1 67000 

16 ■xono 

1 64000 
166000 
1 6 6 000 
16 7000 
1 6 8000 
169000 
1 60000 
16)000 

I 62000 
1 6 lone 

1 64000 
166000 
1 66000 

I I 

1 68000 
169000 
1 70000 
171000 
1 72000 
173000 
1 76000 
176000 
1 76000 
IT610I 
1 77000 
17*000 
1 79000 
180000 
181 000 
182000 
183000 
1 m 4000 
1 85000 
1 86000 
187000 
1 88000 
l 89000 
l< 

191000 

19?000 
193000 
194000 
195000 
1 • 

197000 
198000 
1 99000 
200000 
201000 
202000 
203000 
206000 
205000 
206000 
207000 
20 8000 
209000 
210000 
2 1 1 000 
212000 
213000 
6000 
215000 

2] 60m 

217000 
218000 
2 19 one 
220000 
221000 
H 

223000 

226000 


***• 


««*• 


«*** 
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2 

• . 

•|M| ST/2H24 

. 1 

. 

• • 

1 

?s. 

1 3 

X.6121 






224000 

1 i •< 

Form A l 

. 

H U 

L*- 1 

TM* 

\J 

41 

MA 

TB IX /( 6 

F2(' 

• 

A ) 1 



726000 

620 

H iR RAT 

• , • 

6/ 

30/4 

. 










2270O0 

6 SO 

Fi iRf AT 

1 /SX, 20N I 

1 T 

UL 

STmT 

•» 

vF 

r.T 

. 

» 





?? MOOT! 

640 

^ BRAT 

■ . 

IT 

1 al 

CnVAR 

1 AM 

r.F 

1 • , • ■ 

, 6* 

1 

PC04 IR • 

, 1 

|*I1 

7291)00 

2 

) ) 














229100 

6S0 

FORK’ AT 

I //2*>x, 1/ < 

OF 

FOKF 

r.oA 

< 

T// 

» 







7 VT0OO 

6S5 

F 1 1 R R A T 

1 /*) •', 1 7F'T* 

A i“ 

SI T 1 

no~ ' 

A 

TR 1 

»/ 

1 6iv i.nn 






? R 1 OOO 

660 

FORMAT 

1 ///'SX, I lM 

mF 

T -M 

r n*, f, 

T 

//) 








? R700O 

670 

FORMAT 

1 //29X , |OH 

af 

TFM 

MORI.' 

/ 

/ ) 








? x loon 

4 MO 

FORMAT 

( //26X, 1 |M 

60 

MO F 

R U | (R 

V 

//* 








7 >4000 

690 

format 

)// 

mF 

Tr* 

Ml IR r . 


IT 

H 

r.ORM FHHOR 

s/ / 

) 




714000 

700 

format 

( />x , 2400 ■ 

• 1 

T At 

FLF I 

F 

MTS 

M 

II I lAL/M 1 

Rx. 

A 

h ) /6<V1 

. 1 

1 I 

7 R 4000 

/in 

FORM rtT 

I //»x« 

X , 

6MI- F 

Am * 

, 

F2r. 

• 

. , 

» , 

F 


,6 

Ho P *,07/ 

7 R 7000 

1 


Ml' 

J FF 

i'CFN 

T 

A r.F 

F 

KFOOFMf.Y// 

4 X , 

• 

» |MUS 

’ 

F. '. 

2 3 1 0O(i 


2 1P2E12.4) 

f/o him at ( 4x. 1P3F 12.4 ) 

.1 i K»lPtl2*4»' PUIS |MP ( 
74u Format l/JI 13X, AS l/jnjl .1)) 

Eo»0 


.-.i 


749000 

240000 
? 4 1000 
242000 
?m >ono 


SUBROUTINE" MT»*H(*,P,t,EVl 


IMPLICIT REAL«8(A--I ,□-£) 

OUttHSION Xt 6) .P (6.61 .£ (6) ,EV46 >6) • PX 4 6,3) 

DIMENSION 8*6.6) 

CA* L CUNVET (P.RiXlil.O.PXl 
03 90 I ■ 1.0 
03 40 J - 1.0 
cvt ;. ji ■ pci • j) 
ail.J) *F X 4 I . J) /• 000304 04 42 
40 C3NTIN,'E 
50 CONTINUE 

C A _ L ElbENICV.E *6.1 I 

■ RITE 4 6.600) Ifttl). 1-1.6) 

WRITE 16.610) TEPXII . J) . J*I .61 .!•! .6) 

■ RITE 46.640) 44941 .3) • J* 1 .6) . I ■ l .6) 

WRITE 46.620) UP < I . J) . J -1 .6 ) .1 -1 .6 )? 

■ RITE 4 6. 6 JO ) 4E4I) «I >1 .6) 

RETURN 

600 F3RMAT 4 /5X. 1 2H STATE VE C TOR / ( 602 1 • 1 3 > ) 

6 to ptrinat i nsx . sswcdv a pt aince mmruocut tangent i/iou 2 

620 FORMAT ( / 5X « 2 8HCQ WAR I ANCE MATRIX 4 I NERT I AL >/ 4 602 1 • 1 3 ) 
630 F3RMAT 4 SX . 1 IKE I GEN VALUE S/46071 • I 3) ) 

640 DRMAT ( /5X » 3 JMCO VAR I ANCE MATRIX 4LOC TAN — FD/4602 


1*13)1 


FUNCTION NSAMP4M.0Et.TA) 


C M S.»TS CONFIDENCE LEVEL 

C IF M*IAEVEL I? *9 

C IF M-2. LE VEL IS .95 

C IF M«3. LEVEL I S *96 

C IF M-4» LE VC L IS .99 

C IF M-S. LEVEL I S* 99S 

C IF M«6. LEVEL T S* 999 

C DE_TA SETS UNCERTAINTY 

IF <M-t) 3.2.3 

2 A * 1 * 202 
03 TO 13 

3 IF I 4-2 ) 5.4.5 

O A»t ,049 

03 TO 13 

9 IF4M-3I 7.6.7 

6 A-2,094 
GO TO 13 

7 IF IM-4) 9.6.9 

« ftwtVS T*~ - 

GO TO 13 

9 IF4M-9) 11. 10. II 

10 A-2 • 9 ~i* 

03 T3 |S. 

11 A-3*4)V 

13 N SA WP»t A /t 2 . *0 eLTA >TT» t A /T2* 90EX.TX7 1 
RETURN 
END 


30 


c jcoooo 

REAL FUNCTION U ARN1 • 8 ( I , I K£ V • I F H N , SO ) JOiCCO 

c jC20co 

I** *L I C ! T JO’OOO 

C — - — — — — JO 4000 

C S3-- — — - Tnc DC SI WE 3 ST A '4 JAM J OEVIaTICN 305000 

C *KtAN--“ THI OE SIREd ME A N 306000 

C H THE IX)Ojl*T| JK SIZE JC 7000 

UATA AMLAN/COO/ 30*000 

data iHert/iiTer/ 309000 

DATA H/JCaOf/ JlOOfC 

IE ( |KCV)Si(i4 31 icon 

A I ME *E* IEf<N 3 I 2 0 C 0 

5 1FCI)6. 7 *7 313000 

6 CALL 6 A L SI I mlWI , SO * A ML AN , VAL ,N I 314000 

IF<N»IMtFf. 315000 

03 TO 8 316000 

7 CA.L RANUUC iMtRE . IF RN.VALI 317000 

|HE(iL*If RN jihODO 

B UA RN l ■ VA l 319000 

RETURN 320000 

321000 


c 322000 

SU9R0UT1NE OAUS S( /I »/ *3 .A M, v,M» 323000 

c 324000 

IMPLICIT MEAL* d C A— H *0-2 I 32SC00 

** H 326000 

A«O t OUO 327000 

03 50 I*1,K 320000 

CALL RANOUC I X. I V, VI 329000 

|X«IY 3300C0 

50 A»4*V 331000 

M3*H/12. J J 20 CO 

H2*M/2« 333000 

V«C S«CA-H2) 1/DSQRTCHUlAAM 3J4000 

4ETURN 33SC00 

6*0 336000 

C 337000 

SUBROUTINE RANJU/ll/,1 Y,Y*Ll 33BC00 

c 339000 

IMPLICIT REAL* dCA-rl ,0-2) 34COOO 

DATA JJJ5/1027/ 341000 

IT * I X* 3 J J 5 342000 

IF ( I Y I 5*6,6 34J000 

5 lV*IY+21474e36*7M 344000 

6 VFl. * I Y 345000 

VFL *VFL* *465661 33-9 346000 

RETURN 347000 

E40 348000 


C 3*9000 

SUBROUTINE t I Gt N ( AA , VALU, NR • 4) JjOCOO 

C 351000 

IMPLICIT BEAL* 8 (A-H,a-2) 352000 

REAl*B IND 353000 

C EIGENVALUES AND EIGENVECTORS OF A REAL SYMMETRIC MATRIX 354000 

c 355000 

c 356000 

DHtNSION A C 8*B) ,8(8,8) , VALU(B) ,01 AGIO) .SUPLRul 7 ) , DC 7 ) , V ALL ( b I 357000 

I, S( 7I.CC 7) ,0 ( b) .1 ND C 9) , UC 81 .DUMMY C94) , A A C 64 ) 354000 

E3 U l VALENCE C DI AG C 1 I , J UMM V C 1 ) ) • C SuPERD C 1 ) , DUMMY C 9 ) ) . 359000 

1 C VALLC 1 ) ,D C I ) .DUMMY C 161) , CD Cl ) , SCI ) .DUMMY C 24 ) ) . C bC 1 . 1 ). OUMMYC 31 ) ), 360000 

2 C INUC 1 I ,UC 1 ) ) . Cl I .MATCH) . C TAU.dETAI , CP. PRODS ) • CT .SMALL J). 361000 

3 C A NORM , A NORM 2 ) 362000 

S3R TC X ) «0 SUR T C X I 363000 

SINC V IsDSINC VI 3E4000 

C3SC 2 ) *OCOS( 2) 365C00 

AbSC A ) =0A J SC A ) 366000 

c 367000 

C CALCULATE NORM UF MATRIX 366000 

c 369000 

M«MR 370000 

3R4A * C. DO 371000 

•> *» 372000 

03 1 Ul .N 373000 
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non n n o one* non 


3 R M A a QRUA «A A I I 

1 J»J W*l 

03 2 I ■ I • N 
NI*N*( 1-11 
03 2 JtliN 

1 J “N I ♦ J 

2 t|j||| I A A ( I J ) /OHM A 
J A'O^M^af . njO 

4 03 6 ! »| (N 

5 i)J ft J*1»N 

ft AN JRM 2-ANORM2*A (I ,J)*«2 
7 AN JR** ■ SUP T (AN3HM21 

Ofc.Nfc.HATk lot ST I TV MATRIX 

9 IF (M) 1C, 45, 10 

19 03 40 1 a | ,N 
12 03 40 J*1,N 
2C |F( I-J) 35, 29. 33 
23 H( I • J t a 1 ,0O0 

39 03 TO 4C 
33 U< 1 , J taC.AOO 

40 CONTINUE 

PERFORM ROTATIONS TO REDUCE MATRIX TC JACOBI FORM 

43 IE X I T a I 
30 NN»N-2 

32 IF (NN» 690, 170, S3 

33 03 160 I a| | NN 

60 I I ■ I 4 2 
65 03 160 Ja| I ,N 

70 Tl aA| I , !♦!) 

73 T 2 * A ( I, J) 

OC G3 TO 900 
90 00 103 A ■ I • N 

93 T 2 aC US* AIK, I4|| 4 SUN 4 A ( K , J 1 
IOC At <, J |aCUS*A( K, JI-3UN4A (K ,1 4 1 I 
105 AtX, 141 ) *T2 
110 00 123 K*1 • N 

115 T2-CUS4AI I4| ,K)4SUN4A( J,K) 

120 AU.Kl *C US* A tJ,K| — SUN* A (14 1 ,M 
5 At 1 41, K )aT2 
1 IF (Ml 130, 160, 130 
li - 03 150 K * 1 • N 

135 T 2*CU S* t» ( K , 141)4 S0N4U ( K , J I 
140 at K . J 1 aCUS4B( K , JI-SUN4U (K ,1 41 ) 

150 0( < • I 4 1 ) * T2 
166 CONTINUE 

4JVE JACOB I FOR ■ ELEMENTS A hO INITIALIZE EIGENVALUE BOUNDS 

170 00 200 1*1. N 
180 OlAGt I)«A( |,I ) 

190 VALUt I I aANORM 
200 VA.Lt I 1 a-ANORM 
210 03 230 1*2, N 
220 SU>ERO( I-1).*A( 1-1 ,1 » 

230 3( 1-1 )a( SUPERD(I>n.)4*2 

DETERMINE SIGNS OF PRINCIPAL MINORS 

235 TAUaOaOOO 
240 1*1 
260 MA . C H = 0 
270 T 2 a 0 a OO C 
275 T 1 a 1 a OO 0 
277 03 450 J*1,N 

280 PaOIAGt J J-TAU 
290 IF ( T2 1 3C0 , 330, 300 
300 IF(Tl) 310, 37C , 310 
310 T*P4TI-Q( J- 1 )*T2 
320 G3 TO 4 If 
330 IF(Tl) 335, 350, 350 
335 TIa-1.000 
340 T*-P 
345 GO TO 410 
350 T 1*1.00 0 
355 T*» 


374000 
J 15000 
3780C0 
377000 
378000 
J 7 9000 
380000 
381000 
362000 
362000 
384000 
383000 
386000 
367000 
366000 
369000 
390000 
391000 
392000 
393000 
394000 
395000 
396000 
397000 
3980C0 
399000 
400000 
401000 
402000 
403000 
404000 
403000 
406000 
407000 
408000 
409000 
410000 
411000 
412000 
412000 
414000 
415000 
416090 
417000 
416000 
419000 
420000 
421000 
422000 
42J0f0 
424000 
425000 
426000 
427000 
428000 
4 2 9C00 
4 30000 
421000 
432000 
433000 
434C00 
435000 
436000 
437000 
436000 
439000 
440000 
441000 
442000 
443000 
444000 
445000 
446000 
447000 
448000 
449000 
450000 
451C00 
452000 
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nnn 


J'C GJ TO 4 1C 45J00C 

370 If ( J- I)) 3HC, 350. 380 454000 

380 IF I T2 I 40 8, 39C. 390 455000 

3^0 T *- I • CO C 4SC0C0 

395 GJ TO 4 1C 45 7000 

4 ffl T * 1 ,000 4 1 fcOOO 

C 459000 

C C 3 UN T A^TLMLNTS IN SIGN 480000 

C 461000 

410 IF(Tl) 425. 420. 420 462000 

420 IF I T ) 440, 43C. 430 463000 

425 IF < T 1 4 20, 440, 440 464000 

430 44 TCHSMA TCH4 1 465C00 

440 T 2« T | 466000 

450 T 1 * T 467000 


C fc ST A tiC I SM T | GHTEP HOUNDS ON E 1 Gt N VALUES 

C 

460 03 530 K>1,N 

465 IF | <*MATCM ) 470. 4 7C , 52 0 

470 IF I TAU- VALUlKM 5 30 , 530 , 480 

480 VA.LlXMTAU 

490 G3 TO 53f 

520 IF( TAl>- VfLUIKI) 525, 530, 530 
525 VAlU(K)<TAU 
530 CONTINUE 

540 IF I VAL U< I l-VALLl I ».-5« OU- I 61 570,570.550 
550 IFIVAIUIIII 560. 580. 560 

560 IT I A JS( VALL 1 I 1 /VALUII -I «0l»0) -5*00-16 1 570 .570,580 
570 |>I»1 

575 IF! I-NI 540, 540. 5 90 
500 TAU*< VALLI I IfVALUd 11/2.000 
505 GO TO 26P 
C 

C JACOBI EIGENVECTORS U» POTATIUMAL T R I ANG J L Aw 1 1 AT 1 UN 

C 

590 IF (M| *93. 890, 593 

593 IE X I T *2 

595 03 610 l * 1 , N 

600 03 610 J*l,N 

610 All . J nC.PDf 

ei5 03 E 50 1*1, N 

620 IF <1-11 625. 625. 621 

621 IF < VAL Cl I- I I-VALUI I 1-5. 00-1 4 I 730.730 .622 

622 IF ( VAL U< I — 1 I I 623. 625. 62 3 

623 IF ( A3S( VALUI I 1 /VALUII-1 I -I • OOO) -5*CO-I4 I 730.730.625 
625 CUS*1*00C 

(20 sus=c.oor 

630 03 7C P J * 1 , N 
635 IF(J-l) 680, 680. 640 
640 G3 iO 90C 
650 SI J - 1 ) *> SUN 
660 C( J- 1 I *CUS 
67C 0( J- I .*«T14CUS4T24 SUN 

680 Tl*(OIAG( J I- V*LU( I 1 I AC US- BE TA4 SUN 
690 T2 * SUPtRD I J I 
700 BE TA sSUf-ERO I Jl*US 
710 0(N)*T| 

720 03 725 J«l,N 
725 iND(j)*r*oor 
730 SMALLO * A NORM 
735 D3 780 J*1.H 

740 IF! I0INTI1 NO I J ) 1-11 750,780,780 

750 IF <AuS I SMALLO 1-ABS (00)1)760. 780. 76C 
760 SMALLD-DU) 

770 NN*J 
780 CONTINUE 
790 INO(NN) *1*000 
800 R43DS*l*CO0 
805 IF I NN— I I 810, eSG. 810 
810 03 840 K *2 • NN 

820 1 1 *NN ♦ l — K 
830 41 I 141. 1 l=C I I I l*FR30S 
840 R43DS--PROOS4SI I I I 
850 A( 1 , I ) *PROD S 

F34M MATRIX R93DUCT OF ROTATION MATRIX KITH JACUfal VECTOR MATRIX 
8 55 03 e85 J * 1 , N 


469000 
4 7C000 
4 7 1 00 C 
472008 
4 73000 
474000 
475000 
476000 
4 77COC 
478000 
479000 
480000 
48 1000 
482000 
483000 
484000 
485000 
486000 
48 700P 
488000 
4M9C00 
49C0P0 
491000 
492000 
4.20C0 
494000 
495000 
496000 
497000 
498000 
499000 
5C000Q 

50 1C 00 
502000 
50301.0 
504 0C0 
505000 
5060 CO 
507000 
508000 
5090CC 
510000 

51 1000 
512000 
513000 
514000 
515000 
516000 
517000 

sieoon 

519000 
520000 
52 1000 
522P00 
523000 
524000 
525000 
526000 
527C00 
528000 
529000 
53OC0O 
521000 
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860 03 665 K a I • N 
86 8 

870 03 e88l I a | ,N 
8 75 A( | , J ) a C • CO C 
880 03 ee 12 K *1 iN 

A ( I , J I «b ( liKMUlKItAII • J) 

8852 C3NTINUE 
8851 C3N T INUt 
885 C3NTINUE 

03 net 1*1, N 

Nt»N4( I- I I 
03 886 Ja| ( N 
IJaNlAJ 

886 AA( 13 |aA ( J, I I 
89C C3 X T I NUE 

03 891 I ■ 1 « N 

891 VA.UIII a VALUllltOHMA 
RETURN 

CALCULATE S|N£ AMO COSINE Of- ANGLE 

900 |F (T2I 910# 540* 910 
910 TaSQRT ( T1**2»T2*«2 I 
920 CU3*T 1 /T 
9*. 5 SUN »T2/T 

930 GO TO (9C#6SC), It*IT 
940 GO TO (160,5101. IEXIT 
RETAIN 
EMO 


532000 
5 33000 
534000 
5 JSCCO 
S3C000 
537000 
538000 
539000 
540000 
541000 
542000 
543000 
544000 
545000 
540000 
54 7000 
548000 
549000 
550000 

OF MUTATION 55 1000 

552000 

553000 

554000 

S5S000 

555000 

557000 

550000 

569000 

560000 




C 

SUSROUT 1 Nt BURN ST ( X TR , VTR ,XT • VT • W2 .082 ,WC2 , ALPH.TV • OT, TMRUS. YA. G*<U 
2 > 

C 

C This subroutine propagates the STATE VECTOR THRU The tXJRN 
C WGT IS INPUT IN LBS AND WILL BE 01 VIDEO BY G 

C THRUST MUST tiz IN SAME UNITS AS WGT 
IMPLICIT RE AL* 8 ( A— N • O- 2 1 

DIMENSION XTR(3I#VTR(3I • X T ( J ) , V T ( 3 > • T V ( 3 I , HT H(3I,CTR(3I 
CA.L CROSS! XTR , VTR ,HTR 1 
CALL CKCjSS(MTR.VTR,CTR> 

CM A G a FNORM(CTR) 

VMAi* a F N 1PM ( VTR) 

MM A G ■ FNORH(HTR) 

10 Alpha a ALPH a.Cl 74S3293D0 
BETA * 1#57C 79 6 3 2 73 0 4 ALPHA 

C0SAL a OCO S ( ALPHA ). 

C3 SUE a DCOSIBETA) 

YAW aVA a. 01 745329300 

SYAW a DSINIYAWI 
CYAN a OCOSIYAwl 

C CONVERT LBS TO ( KG-KM / SEC»«2I 

Ml a w2 4 #00444822 1 7DC 
DWl «OW24 #00444822 1 730 

tC|a«C2 4 #00444822 1730 
THRUST ■ TMRUS 4 • 0 0444 82 2 1 7D0 

03 I 1 a 1,3 

1 TV! I I 4 THRLST4 (C0SAL4C YAW4VTRII l/VMAG 4 CCS BE4CY AW 4CT R( I 1/ CM AG 

1 4 SYAWAMTR ( I l/MMAG ) 

G a 9# 806650- 3 
W a wl/G 
OW a OW1/G 
WC a WC1/(G43T) 

TEMP a OL3G! 1#004 WC 40T/ ( W43 Ml) 

R MAG a FNORM ( XTR I 
CONTINUE 

03 2 1 > 1,3 

XT! II a XTR ( | l-GMLf* XTRII ) 4I0T442 |/(2.O0 4RMAG443) 4VT R( 1 I40T4TV! 1 14 
1 ( ( W40W4WC43TI/IRC442) 4TEMP-3T/WC) 

VT (II a -GML4XTR! I | 40T/ (RMAG443 14VTR! I |4TV( I I/WC4TEMP 
2 CONTINUE 
RETURN 
END 


56 1000 
562000 
563000 
554000 
565000 
566000 
567000 
568000 
56 9000 
570000 
571000 
572000 
573000 
574000 
575000 
576000 
577000 
578000 
579000 
580000 
58 1000 
582000 
583000 
584000 
585000 
586000 
587000 
588000 
589000 
590000 
591000 
592000 
593000 
594000 
595000 
5960C0 
597000 
598000 
599000 
600000 
601000 
602000 
603000 



C 604000 

SUBROUTINE PARTAL ! PAR , TR , TV# T50 • ALPH ,BET • T MRUS I 605000 

C 60600C 

IMPLICIT REAL48IA-H.O-ZI 6C7000 


34 


UJJULt PRECISION PAPO.3). TSO (6.3) 
DIMENSION THIS) » TV ( SI • Vlll. /(J) 

T ■ THMLSA*F044*a22| 700 
A_>HA ■ ALPH /S7.29S779SI JO 
UtTA • btT /S7*2 VS77VSI 30 

iA.PHA ■ OSIN (ALPHA) 

CA.PHA ■ OCO SI ALPHA ) 

Sat TA • CSINIDCTA) 

CbtTA * OC OSIULTA) 

•*4(1.1) ■ CALPHAA cbeta 
PA- 4 1 2.1) • - SALPHA'Cbt TA 

UK 3.1) • SttETA 

PA4(I.2> ■ - T*5ALPHA*CMETA 

3 A 1 1 2. 2) • -TAP AH (I .1 ) 

PAR (3.2) a 0*00 

PARI 1.3) ■ -T*L ALPHA* SdETA 

PA4I2.3) ■ TA SAlPHAA S dt TA 

>1)13.3) a TACUETA 

VM ■ FNORMI TV) 

CA.L CROSS I TR .TV. 2) 

CALL CROSS (2. TV. T ) 

TM m FNORMI V ) 

2M ■ FNOHMI2I 
03 I 1 a | . 3 
TSOI I . I I > TVI I l/VM 
TSOI I . 2 ) ■ VI I ) /TM 
TSOI 1.3 1 • 2(1 ) /2M 
1 CONTINUE 

AET1PN 
EMD 


6 CeOOO 
009000 
O I 0000 
A) 1000 
E12000 
013000 
<14000 
OISOOO 
41(000 
61 7000 
O 1 0000 
039000 

021000 
023000 
02 3000 
024000 
028000 
02O000 
027000 
02*000 
029000 
030000 
031000 
032000 
033000 
t 34000 
039000 
039000 
040000 


SUBROUTINE POVERXIXT. V T • XL T . XMT • XMOOTT ,T T • GMBOU .TSP.P.PMl 


THIS SLUHGLTINL COMPUTES THE STATE ThANSITION MATRICES P.O.R 
AT TIME T 3 UR 1 NO PUMEREO FLIGHT. AND THE FOPMARO PROPAGATION 
MATRIX PM* 


IMPLICIT PL ALAOIA-H ,A-2> 

DIMENSION XT! 3) *VT( 31 • XlT I 3 ) • XTR I 3 ) .VTRI3).LTP(3).P(6.6).PM(6.3) 
03 30 I a 1.3 


XM(I) « 
VT4I1) i 
L T 1 1 I ) i 
VT4 aXMT 
M03TTR i 
TR a TT 
gmbooy ) 

T A TSP 
RSQ « I 
RCU8E a 


XT! I ) 
VTI I I 
XL TII ) 


X TP I 1 I A A 2 ♦ XTR ( 21**2 ♦ XTR(3)**2 I 

RSQ A OSORT(RSO) 


OT a T - TR 
A A GMBODV / ROUE 
8 a A A OT 
C ■ B * UT/2.0D00 
0 > 3* 00 CO AC/ RSQ 

E > OT / RCUUE 
P A t * DT/ 2*0000 
G > 3* OU CO • U / HSU 

H a M TR ♦ MOOTTR A DT 
* ■ DLOG1H/ MTR)/ MOUTTR 
03 1 I > 1 • 3 

03 2 J ■ 1.3 

*11*31 ad* XTP(l) * XTR I J) 

P( 1 .3 *3 ) > 0*0000 

PI 1*3.3 > « G * XTR ( I ) A XTR I 3) 

P( I O.JOI « 0* 0000 

IF I I *NE • J ) GO TO 2 

P(I.J) * Ptl.J) < 1*00C0 - C 

PI 1.3*2) « OT 

PI 1*3.3 ) * P(I*3.3) - 8 

PI I*3.3*3) a 1.0000 

C3NTINOE 

CONTINUE 

MAG a OSORTILTRIl ) * *2 *L TR 1 2 ) A *2 *.LTR I 3 > *#2 ) 
<1 A MDOTTRAK 

XTERM a -IOT-(MTR/MOOTTP*OT)AK1 )/MOOTTR 
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X TERM I ■ ( 2.0300* 3 T- ( 2 • ODO 0* MTR/MOOT Tft*OT I *K I ) / C MOOT T R *M AQ* *2 ) ttaCOO 

03 in 1*1,3 ewoooe 

03 10 J a 1 , 3 tvOOCO 

»«l I ( J » aXTERM|*LTR( I t • L TH ( Jl 691000 

IMI.rQ.J) P M ( I • J ) • PM (I ■ j)« xTC fix 642000 

|0 CONTINUE 692000 

XTE4MI • -(* t/MAG-MDUTTR*DT/(MAG*H) )/ ( MOOT T R*maG) 644000 

03 20 I >4.6 646000 

03 20 J * 1 • 3 696000 

»W( 1*21 • XT£RMl*LTR( 1-3) *LTR( J) 69)000 

IM 1-JitUtJ) PM 1 1 « J 1 > PM(|,j|»k 69*000 

20 C 34 T l Nut 699000 

RETURN 700000 

LNJ 701000 


c 702000 

SUBROUTINE MTRPL T (A.B.C.NRA, NLA, NC B , NA, Nb « NL I 703000 

C 704000 

IMPLICIT RE AL*P (A-M ,0-2) 705000 

03 UBlE PRL'ISION A(NA,1), i) ( No • I ) • C ( NC , 1 ) 706000 

03 1 I ■ 1 « Nk « 70 7000 

03 I J«I.NCb 7COOOO 

C < I *• J ) « 0,0** 709000 

03 1 X • l , NC A 710000 

C(I,J) * C 1 1 * 31 ♦ All «K ) *B ( X , J) 711000 

1 CONTlNUt 712000 

Rt T URN 713000 

640 714000 


C 

SUBROUTINE 7 DO 00 V( SO ,TAU, MU , PSI ,S,P*PI • PMU . PO MJ , ACC, ACCO . R. RO ) 

c 

c 

C GENERAL SULLTIUN OF T 60 BODY PMlIbLE » ■! TH PARTIAL DERIVATIVES 
C F3RTRA4 4 DOUBLE PRECISION SUBROUTINE FOR I »M 7094 WITH IbSVS SYSTEM 
C SEE APRIL 1965 ASTRONOMICAL JOURNAL FOR FORMULATION bV M, GOODYEAR 
C 

C CALLING SEUUENCE IS AS FOLLOWS 

C CA_L TPOUDYl SO ,TA U , MU • P SI ,S,P,PI , PMU . PC Mu . ACC , ACCO ,R,RO ) 

C 

C DOUBLE precision tlUAN TI ties in calling SEQUENCE are as follows 
IMPLICIT REAL48IA-4 ,0-2) 

DOUBLE PRECISION SO 1 6 ) , TA U , MU , PS I 
1, SI 6)«P(6<6) » 3 I (6,6) .PMU (6) ,P0MU<6) , ACCO) • ACCO 1 3 ) , R, RO 

C 

C INPUTS 

C SP( 1 ). SOI 2) • SOI 3) «X(J ,Y0 .20*P0SI TI ON COMPONENTS AT REFERENCE TIME TO 
C SCI 4) , SOI 5) , SOI 6 ) >XOC.YD0.2D0 «VELOCI TV COMPONENTS AT REFERENCE TIME TO 
C TAU-TIME INTERVAL (T-TO) FROM REFERENCE TIME TO TO SOLUTION TIME T 
C MU*C34STANT IN DIFFERENTIAL EQUATIONS 1 XDD , YOU , 2DD ) a -MU * 1 X • Y , 2 )/( R44 3 ) 
C PSI * APPRO X I M ATI ON FOP FINAL SOLUTION PSI OF KEPLER* S EQUATION 
C 

C OUTPUTS 

C PS I = GENERAL I 2ED ECCENTRIC A NOMA L V* SOLUT I ON CF KEPLERS EQUATION 
C SI 1 ).S( 2).S( 3)*X»V»2 *POSI TI ON COMPONENTS AT SOLUTION TIME T*TOFTAU 
C 31 4 ) • SI 5), SI 6) *XD. YD , 2D «VELOCI TV COMPONENTS AT SOLJTlON TIME T*TO*TAU 
C P( I, J IMPARTIAL DERIVATIVE 3S(I)/DS0(J) OF SCI) WITH RESPECT TO SO(J) 

C P I ( I , J ) aPAR T IAL 3 SO 1 I ) /O S ( J ) WITH ROLES OF TO AND T REVERSED 
C PMUC I IMPARTIAL DS(I)/OMU OF S1I) WITH RESPECT TO MJ 
C PCMUl I ) ^PARTIAL D SO 1 I ) /O MU WITH ROLES OF TC AND T REVERSED 
C ACC1 I )«-MU» S( I )/(RM3) *ACCELERATI ON COMPONENT AT SOLUTION TIME T 
C ACCO! I ) *-MU» SCI I )/(R0**3) aACCELERATI CN COMPONENT AT REFERENCE TIME TO 
C R*RADIUS AT TIME T * SQ UARE ROOT OF ( * 4*62 *■ V4*£ *.l *«2 ) 

C ROaRADIUS AT TIME T 0 a SQ UARE RCOT OF ( XO • *2* YO *42 *2? **£ ) 

C 

C ADDITIONAL DOUBLE PRECISION QUANTITIES FOR COMPUTATION 

2. SI GO. ALPHA .PSI N.PSI P .A .AP.CO ,C1 ,C2 • C3 • C4 ,C5 X3 , S i , S 2 , S 3, CT AU, OTAUN 
3. OTA UP , U«FM 1 .G.FD.GOMl 
C 

C START OF INITIAL COMPUTATIONS 

C C3MPUTE RADIUS ROaSQUARE ROOT OF 1 XO 4 42 F YO 442 + 20 4 42) 

SI -UMAX H DABS! SO(l).) .DABS (SO (2) ) .DABS (SO (3) ) ) 

S2>< SC( 1 )/Sl )**2*(S0(2) /SI) **2* (SO (3) /SI ) 4*2 
40*2*00 
10 4 aft 0 

R0*(RFS2/R )*• 500 
IF(RO*LT*R) GO TO 10 
R0*R0*S1 

C COMPUTE OTHER PARAMETERS 


715000 
n toco 
717000 
7 1 8000 
719000 
720000 
72 1000 
722000 
722000 
724000 
725C00 
7260C0 
727000 
726000 
729000 
730000 
721000 
732000 
733000 
734000 
735000 
736C00 
737000 
738000 
739000 
740000 
741000 
742000 
743000 
7440C0 
745000 
746000 
747000 
74600C 
749000 
750000 
75 IOC 0 
752000 
7S3000 
754000 
755000 
7S6000 
757000 
756000 
759000 
760000 
761000 
762000 
763000 
764000 
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5 1 00 » SO ( I >» SCI 4 )♦ SO (21* SO 15 I ♦ SO (.11 *S0 (6 » 

A^ >m* * SC ( 4 ) • *2* SO ( 5 ) **2 ♦ S0( t> ) #*2 -2 • L»0 •MO/R0 
C IN 1 T I Ah. 1 21 SI « It S MQU COUNT M TU it RC 

M «0 

C INITIALIZE ML UNO S PSIN ANU PSIP FUM F.»I OR StT PS 1*0 If T AU»C 
|F|TAU) 20,30,40 
20 PSIN*-taD»3e 
PS!P*0,0f 
OTA CM -PS IN 
OTA UP*- TAU 
03 TO 5C 
30 ®SI»0,0C 
03 TO IOC 
40 ®SlN»OaOC 

PSIP-41 *0*36 
0T4 UN— TAU 
OTAUP aP SIP 

C USE A®»ROXlMATION FO» ®S! |F IT IS BET «£IN BOUNDS PSIN AND PSIP 
SO |F(PSIaOT a PSIN a ANO«PS| a LT.PSIP| 00 TU 100 
C TPT NE*TON*« MtTMOD FtJR INITIAL PSI SET EQUAL TO ZERO 
®SI aTAU/PC 

C SET »S I *TAU IF NtaTON'S METHOD FAILS 

IF < PS! a LI a PSlN a OR a PSl a 0E a PSI P) PSI > TAU 
C LNO 3F INITIAL COMPUTATIONS 

C 

C BEGINNING OF LOOP FOP SOLVING KEPLER* S EQUATION 
C BEGINNING OF SERIES SUMMATION 

C COMPUTE ARGUMENT A IN REDUCED SERIES OBTAINED BY FACTORING JUT PSI'S 
ICO A aAL PMA 4PSI-PSI 

IF ( OAtiS ( A » a LE a 1 ,D0> GO TO 120 

C SAVE A IN AP ANO MOO A IF IT EXCEEDS UNITY IN MAGNITUDE 
AP ■ A 

1 1C M *N ♦ 1 

A a A 4 # 2 SO C 

IF I OABSI A ) a GT a I .001 GO TO 110 
C SUM SEMES C5X3»3*S5/PSI4*3 ANO C4*S4/PSI**4 

I 20 C SX Ja( | ,of H l a 30+(l a un» ( I » JO* (I a DO * ( 1 a OO ♦ ( 1 , O0 ♦ A/ J4 2 a 00 » * A/ 2 72 ,oO I 
l**/21CaDO)*A/13 6.DOI 4A/I 1 0, DO I *A/7 2 a D0 I 4A/42 a O0 1/40 aOO 

C 4 «( 1 aDCM l, Jft* ( I aD OF (I ,U0» (1 a OO ♦ I l a DO ♦ I I a JO 4 A/ 30 b a DO )*A/240aD0 I 

|4A/|e2aDr)-A/132aO0) • A/90a 00 > 4 A/ 56 a DO I -A/ 30 a DO 1/24.00 

C COMPUTt SER IL S C3*S3/PSI *43 .C2-S2/PSI *42 .Cl -S1/PS1 .CC-SO 
C3-I a SUCTA4C SX3) /3a 0« 

C 2* a SOC 4A4C4 
Cl- 1 aO C ♦ A 4 C 3 
CO* I »D C ♦ A 4L 2 
IF(M.LE.C) GO TO 140 

C 06 4 30 SERIES CC AND Cl IF NECESSARY »ITH DOUBLE ANGLE FORMULAS 
130 C1-C14C0 

C0»2aD0 4CC4CO-IaD0 
K*M- I 

IF(M.GT.C) GO TO 130 

C DETERMINE C2,C3,C4,CSX3 FROM CO,CI.AF IF DEMOO REQJIREO 
C2-ICP- » ?00 l/AP 
C3»«C I - l,b •*» /AP 
C4«(C2-,.un»/AP 
C SX3-I 3aDC4C3-aSOO) /AP 
C C3M ■* UTE SERIES S1,S2,S3 FROM Cl ,C2 ,C3 
140 SI-CI4PSI 

52- C24PSI4PSI 

53- C34P SI4P SI4P SI 

C END 3F SERIES SLMMATION 

C COMPUTE RESIDUAL OTAU ANO SLOPE R FOR KEPLER'S EQUATION 
G-R 04S14SIGC4 S2 
DTAD-I G4MD4 J 3I-TAU 
R -OABSI R04CC4I SI GO* S I ♦ M u* S2 ) I 
IF(OTAU) 2CC.30C.210 
C RESET BOUND 
200 PSIN-PSI 

OTA CM -O TAU 
G3 TO 2 20 
2 1 C PSIP-PSI 

OTAUP -OTAU 

C TRV NEKTON* S ML THOO ANO INITIALIZE SELECTOR N 
22C PSIaPSI-OTAU/R 
N*0 

C ACCE*T PSI IF IT IS BETmEEN BOUNDS PSIN AND PSIP 
230 IF< PSlaGT.PSlN.AND.PSta LTaPSIPI. GO TO I 00 
C SELECT ALTERNATE METHOD OF COMPUTING PSI OR STOP ITERATIONS 
N *N ♦ 1 


765000 
766000 
707000 
700000 
769000 
770000 
771000 
712000 
7730C0 
774000 
719000 
770000 
777000 
770000 
7 7 VO CO 
760000 
70 1000 
702000 
703000 
764000 
765000 
7C6000 
767000 
700000 
769000 
790000 
791000 
79200C 
793000 
794000 
795000 
796000 
797000 
790C00 
799000 
6C0000 
BO 1000 
80Z000 
603000 
B04000 

eosooo 
006000 
607000 
600000 
609000 
010000 
61 1000 
612000 
613000 
•14000 
015000 
016000 
617000 
016000 
019000 
620000 
621000 
022000 
623000 
624000 
625000 
626000 
627000 
620000 
629000 
630000 
631000 
6J2000 
633000 
634000 
63S000 
83C000 
637000 

836000 

639000 

640000 

841000 

■42000 

643000 
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S3 TO I I .2.3.4,300 .ft. 

C T6V UCAtMt NT | NG JO UNO ll T1 0 TA U MANtST 2tRU UT IHI MAT lO 4ALT4U/T4U 
t IF! 0A9 9(DT4UN),LT. jabsotaupi » PS|*PSIN*<I.U0-»A.U0*'>T AUN )/T Au ) 

|F(OAJi(uT*UM.LT,OAa,|UTAUM I PS 1 ■ PS 1 P A ( | • UO - 1 4 . oO AUT 4uP|/TAU) 

63 TO 2 JC 

C thy UJuBlINC. antNU CLOSEST To 2ERO 

2 IFC TAU.GT.O.OP) P SI »PSI N4 PS! N 
in TAU.L T.O.DO > P SI *RSI PA PSl P 
63 TO 2 30 

C TRY IS TIM POL A T 1 ON UtTttLN UUUNUS 

3 PS I *PSINAlPS|P-PSINI A(-DTAUN/IOTAUP-DTAUN| | 

63 TO 230 

C THY HALVING UET6CEN BOUNDS 

A PS I SIN ♦ I P 5 IP-PSI N > •• 900 
63 TO 230 

C tNU OF LOOP ft OH SOLVING KLPLER’S EQUATION 
C 

C COMPUTE REMAINING TMRlE OF ftOUR FUNCTIONS F Ml • 6 • F U. GUM I 
300 FM| •-■U*S?/P0 
F 0 •—MU# S I/R0/R 
60M I «-MU« S2/R 

C COMPUTE COORDINATES AT SOLUTICN T I ML T-TOATAU 
OO 310 I >1 . 3 

S( I I -SO I I IMF Ml* SOI I IaGASOII A3) I 

sc i»?i*(n)*somftsointsomiii»so(ift]i 

C COMPUTE ACCELERATIONS 

ACC I ! )— MUASIt I /R/P/R 
310 ACCOI | »• -MU* SOI l l/RO/RO/RO 
C END 3F COMPUTATION FOR COORDINATES ANO ACCELERATIONS 
C 

C COMPUTATION OF PARTIAL DERI vATIVES 
C COMPUTE COEFFICIENTS FOR STATE PARTI ALS 

U> S2ATAU4MUA (C4-C5X3) APSI APSI APSI APS l APS I 
*1 1. 1 I — IFDASIAFMl/ROI/RO 
P( 1.21 »-F D* S2 

P(2.1)> FM1ASI/R0 

PI 2. HI FM1*S2 
PI I . 3) a PI I .2) 

PI 1.4) — COM I • S2 
PI 2 • 3 ) ■ PI 2.2) 

PI 2.4)* 6* S2 

»l 3. I |— FDA |C O/RO/RA l .DO/R/RM .UO/RO/RO ) 

»l 3.2) — I FDa SI *60 Mi /R)/R 
PI 4. t )»-PI 1 .1 I 
PI 4, 2)— PI 1.2) 

PI 3.3)> PI 3.2) 

P{ 3.4) — COM 14 SI /R 
PI 4. J)»-PI I .2) 

PI 4.4)— PI I .4) 

C COMPUTE COEFFICIENTS FOR MO PARTI ALS 
PI l.SI — Sl/RO/R 
PI2.S)- S2/R0 
PI 3. 9 ) ■ U/RC-S3 
PI I.O— PI I .9) 

*12.6)* S2/R 
PI 3.6) — U/RAS3 
DO 400 I ■ 1 » 3 
C COMPUTE MU PAR T I ALS 

PMUI I I — SI I ) API 2. 3) ♦ SI I 43)4 P| 3. 9) 

PMUl I 43 I ■ SI I )AP( l .3) AS 1 1 43) API 2 .51 AACCI I ) API 3.9) 

POMUI I) —SCtt )4PC2.6)»S0II*3) API3.6) 

POMUI IA3)« SOU )*PI1 ,6)4 SOI I A3)AR(2 ,6) A ACCO I 1 ) API 3.6 ) 

C MATRIX ACCUMULATIONS FOR STATE PARTI ALS 
OO 400 3 ■ I , 4 

Pit J. I )> PI J. II* SOI I I API J.2) AS0UA3) 

400 PIIJ.143)> PI J.3)*S0II I APIJ.4) ASOII A3) 

03 410 1*1.3 
03 420 3*1.3 

PII.J) *SI I I API I I • J) ASH A3) API 12 . J) AU AS I I A3 ) A AC CO I J ) 

>11.343) «SI I )API 1 1 .34 3) ASH A3). API 12 .3431-0 AS I I A3) ASOI J A3) 

PI 143.3) * S I I I API 13.3) ASIIA3)«PI 14.31 AU AACC I I ) A ACCO I 3 ) 

420 PI I *3.3 A3) *SI I ) API | 3. 3* 3) ASH A3 I API 1 4 . 34 3 ) -U A ACC I I ) ASO I 3 A3 ) 
•Tift) APII.I1 AFR1A1.00 
PI I. 1A3 I *PI I . I A3) A6 
PI IA3. I ) API 143, i) AFO 
410 P| lA3,lA3)aPll43.|4 3)A60MIAl a DO 
C TRANSPOSITIONS FOR I N VE R SE STATE PARTI ALS 
OO 430 I ■! . 3 
OO 430 3»1 .3 " " 


644000 

649000 




690000 
691000 
6S2000 
69 3006 
694606 
690000 


fcS/OOO 

696000 

699060 


661000 

662000 

663000 


669000 


067000 
660000 
*69000 
670000 
671000 
672006 
6 1 3000 
674000 
679060 
67*006 
677606 
*16006 
*79000 
6*0000 
661000 
662000 

604000 

*69000 

666000 

667000 

666000 

009000 

690000 

691000 

692000 

693000 


699000 

096000 

697000 

090000 

099000 


901000 

902000 

903000 


909000 
906000 
907000 
906000 
909000 
910000 
91 1000 
912000 
913000 
914000 
919000 
916000 
917000 
910000 
919000 
920000 
921000 
922000 
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s w > u n •* n n r> n » » o> •• n a a r n o n a a 


»li j« pii.jj «2M«0 

MUOitl t-PIIOiJI 92*000 

>1(3.1*31 *-P(t,3*.3) 92 !0C0 

• 30 Pll Jill • »(IMi JMI 92C0C0 

C 190 Of COMMUTATION FOR PARTIAL OlRIvATIVLS 922000 

920000 

t NO or PROGRAM - ACL OUtPOlS H A Vt etc IN CCNPOtfcO 929000 

• I TURN 920000 

t 90 9 J 1000 


SUJROUTIM CONVLT (P.R.V.K.PUI 
K«1 NO TRANSFORMATION 

K* 2 p TRAN 52 TP Mi O KRCM LOCAL TANGLNT TC LuUATUH 90 

K*0 P TRANSFORMED FRCM £ ouatoh so to local tanglnt 

U IS TML TRANS* ORMATI UN 
IMPLICIT SUL#U«A-HiO*/l 

OJUdLt PRECISION P(6,6). 0(6. 61. UUP. Pi.lb.bl 

DIMENSION H(3I. V(3 J • RAV(JI. RXVAH(J) 

ir IK ,tC, II HLTUHN 
CA.L CROSS (R.V.HAV) 

L*.L CROSS (WAV. R, MavaHI 
*9 • r NORM IN) 

H A 99 a r NORM |nl«| 

RAVARN a FNORM(RAVAR) 

03 I la I, J 

3(1.11 • R A V AH III /HAVARN 

0(1.2) ■ R AVI I ) /HAVK 

0(1.3) aRII) /RN 

03 ? J a | , 3 

9(10 .JO) ■ QII.JI 

3(1*3. j I a n s or 

0(1. 3*3) • 0,0 0 

C 39 T I NUt 

C39T1NUC 

IF (K ,C0, 2) GO TO 3 
03 * I > I, 6 
03 5 J« I .1 

0 09 a U| I .J) 

0(1.31 • 0(3.1 ) 

3(3.1) « 0 LM 

C39TINUL 

CALL MTRX (Q. P. PQ . 6. 6.-1) 

RETURN 

£90 


9J2000 
9JJ000 
934000 
9 JS000 
9J60C0 
937000 
93*000 
9 J9000 
9*0000 
9*1000 
942000 
94 J000 
94*000 
949000 
946000 
947000 
94*000 
949000 
9SOOOO 
991000 
952000 
953000 
954000 
955000 
956000 
957000 
95*000 
959000 
960000 
961000 
96 20 00 
963000 
9(4000 
*65000 
966000 
967000 


SUBROUTINE M TP A|A .3 .Q.NP . NC • M ) 

IMPLICIT RC AL*8 (A-H ,0-Z) 

L> J UtlL L PR CC I SION A.J.C.O.U 
01 MENS I UN A(6.6),J(NR.NC> .C (6 .6 ) .0 
*■0 G! VL S QaAU 
Ma-1 GIVES Q «AU A TRANSPOSE 
03 7 I » I .6 

D3 I 3*1. NC 
0( 3 )aP,0C 
03 I K a I .NR 

0(3) *C ( 3)*A( I «K)*d(K,3) 

IF(M) 4.2.4 
CONTINUE 
03 3 3»t.NC 
C( 1.3) *0 ( 3 ) 

G3 TO 7 
03 5 3*1,6 

C( I ,3 ) «CaO0 
00 b K ■ 1 ,NC 

C( 1,3 l*C( 1 .3 )*3(K)4A(3,K) 

C3 N T I NUt 
03 8 1*1.6 
03 8 3*1.6 

8 3( I .3 ) *C( I .3 ) 

RET URN 
END 

C 

SUBROUTINE T Ad 1 ( A . S .NOVAR ,U0 O ,F RE G 


968000 

969000 

970000 

971000 

972000 

6 ) ,0(6 ,61 973000 

914000 

975000 

976000 

977000 

970000 

979000 

9*0C00 

901000 

9*2000 

9*3000 

9*4000 

9*5000 

9*6000 

957000 

9**000 

909000 

990000 

991000 

992000 

993000 

994000 

995000 

996000 

PCT.STATS.NU.NV 1 997000 
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c 



REALM* 2(41 

9N«000 
• 999000 


DIMENSION At 1 )» Sill • UtJli I 1 » .MteOI 1 I 'PCTI I 1 .ST ATI ( 1 ) 

1000000 


DIMENSION iUOIJI 

ICOIOOO 


00 S 1 • 1 * 1 

1 OC 2000 

5 

■ S3 f 1 1 ■ UBOt 1 1 

1003006 


VM IN » uec 75 

1 OQ 4000 


VMA A ■ - I. PE 75 

teosoeo 


I J a NUM INOVAR- | 1 

1006000 


03 30 J ■ 1 'NO 

IC07000 


IJ • IJM1 

ICOOOOO 


IF (SUM t C • 30 .10 

1009000 

10 

IF lA(UI-VMlN) 16.20.20 

IOIOCOO 

15 

VM | N • A 1 1 J 1 

101 1000 

20 

IF <A||J)-VMAX| JC.30.25 

1012000 

25 

VMAX • A 1 | J | 

1C 1 3000 

:0 

CJNTlNUt 

1014000 


STATSI4I a VM | N 

1015000 


STATS! 51 ■ VMAX 

1016006 


IF 1 UBOt 1 1-LBOI 31 1 4C '35 '40 

1017000 

35 

Ot*J( II ■ VM | N 

101 6000 


DO 313 1 a VMAX 

1019000 

40 

|NN ■ UUGI2I 

1030000 


D3 45 I a 1'INN 

IC2 1000 


FREQI II a OaO 

1022000 

45 

PC T 1 II a Oaf 

1023000 


03 50 1 a 1 , J 

1024000 


21 I I • CaDO 

1025000 

50 

STATS! I I * OaO 

1 026000 


SINT a AOSI 1 UHD 1 3 1- UBOt 1 1 1 /IUUUI2I -2a0 1 ) 

1027000 


SC N T a 0.0 

1026000 


IJ • NO • 1 NO VAR — 1 ) 

1029000 


03 75 J > | ,NJ 

10J0000 


IJ ■ IJ*I 

10JI000 


IF tSIJI) 5*. 75, 55 

1 032000 

55 

SC N T a iCNTFlaO 

1 0 J 3006 


STATS! 1 1 a STATS! II FA <1 J I 

10J4000 


ST ATSt 31 a ST.\TSU|FA(I JI*A(I Jl 

1035000 


21 41 a All Jl 

1036000 


2111 a 21 1 IFZI 4 1 

1037000 


2131 a 2i3l*2t4l**2 

1036000 


TEM» a UOOIII-SINT 

1039000 


IN TX a INN-1 

1040000 


03 60 1 a i.InTX 

1041000 


TEMP a TEMP ♦ SI N T 

1042000 


IF <A( IJI-TE MP | 70'6f.60 

104 3000 

60 

CJNTlNUt 

1044000 


IF <AtlJ|-TLMP) 75 '65 .65 

1045000 

65 

FREOflNNI a FRt U 1 1 NN 1 ♦ 1 a 0 

1046000 


63 TO 7 « 

1047000 

70 

FRLDt 1 1 a FREQ! 1 >'l .0 

1046000 

75 

CONTINUE 

1049000 


STATSI 11 a 2111 

1050000 


STATS! J 1 a 2131 

1051000 


03 60 1 ■ 1'INN 

1052000 

00 

PC Till ■ FRfcOll |A lOOaO/SCNT 

1 OS 3000 


IF < SCN T- 1 a C 1 65.65,90 

10S4000 

65 

STATS! i 1 - OaO 

1055000 


STATSI 31 a CaO 

1056000 


03 TO se 

105 7< 50 

90 

ST ATS! 2 1 ■ STATSI 1 1 /SCNT 

1 C 5 6> : a 


STATSI3I * SORT ( Alls II ST AT S! 3 1 - ST AT S ! 1 1 MST ATS i 1 1 FS CNT l/t S CNT - I 

laOIII ICST'O' 

95 

03 100 1 a 1,J 

1061 /Of 

100 

US3I I 1 a WU0I1 1 

ICt lOw-J 


RETURN 

10. 20C' 


ENO 

10 3030 


c 1064000 

«t*L FUNCTION 00T»8 ( X • VI 1065000 

C 1066000 

IMPLICIT RE AL* 6 !A— H *0—2 1 106T000 

DIMENSION X(JM(3I 1066000 

DOT > AC 114 Till ♦ X(2I»V(2I 4 Xt3)AV!3) 1069000 

RETURN 1070000 

END 10JI000 


1072000 


40 



PEAL HjW TIUN (N.A.Tl 

c 

I '••Lie IT «f *L» 8 IA-M , 0 - / I 

C CJMPUTlS 4-UUAJHANT ARCTANv#LNT UF T/A IN RADIANS 

C Nal6« ANGLE LltS IN R ANbt I C . I NO I DM# 

C N • I mO ANQLl LUS IN RANgL t - 1 <*0 • t (• 0 » UL<# 

r»lM*<l)lllMM rssBftjp 

AA aOAiiSI A | 

TA •OAb SI T | 

IP* AA-VAI 1.1.2 
» l a A / TA 
00 TO 3 

2 /■ T/AA 
TA «AA 

3 D «D SQM T I I . 90 / I 

TA «TA*0* A 

IM TA IA.A.9 
A AAA TN S ■ TP I / 2 • OU 0 
GO TO 6 

5 A P* TNS s 2. OO P*DATANI T/TAI 

o imnmkit.t.t 

r IFI ARkTNSIE. 9.9 
8 A A A TNSaARATNS*T*| 

V R t T UN N 
END 


10/3000 
IO/AOOO 
10ZS000 
I 0/6000 
IC//000 
101*000 
I0/V000 
1080000 
10*1000 
10*2000 
1002000 
10**000 
1000000 
10*6000 
I CO 1000 
10*0000 
IC09000 
1CV0000 
I 0 V 1000 
1 0 1*2000 
1003000 
1CVA000 
1009000 
ICVAOOO 
I CO 7000 


C 180*000 

UAL FUNCTION F NORMA BIXl 1099000 

C 1100000 

IMPLICIT PEAL* 8 IA-H .0-21 1101000 

D I TIN SI ON a I 31 1102000 

1 FNDRM aOSOPTlAl | ) *• *♦ A( 2 I ••*» A 1 3 1 4*2 I 1103000 

3 PE TURN 1 104000 

END 1109000 


C 

c 

c 


SUBROUTlNt CPU SSIA.B.CI 

IMPLICIT UALMIA-H.O II 
DIMENSION A I J) .U( 3) «C I 31 
0 I MEN S I ON At 3) *B I 3 I , C ( 3 ) 
Cl 1 l-AI 21*61 3I-A1 J)«tf(2l 
Cl 2 I “A | 3 )** I I 1-A| 11*0(3) 
Cl 3 I -A I 1 1 ABC 2I-AI2I *UI1 I 
RETURN 
END 


I 1 CEO 10 

I 10 7000 
1 10*000 
1 109000 
1 1 1 0000 

I I 1 1000 
1 112000 
1 113000 
1 11*000 
1 1 19000 
1 1 16000 


c 

SUBROUTINE ORO I A *D A . L .OE I 
C 

IMPLICIT REAL*8IA-H ,0-21 
DIMENSION OE 161 
DIMENSION At 31 .DXI3) ,Ut3l 

C THE F 3L L 3 ■ INC STA TE ME NT | SI HAVE BEEN MANUFACTURED BT T ML TRANSLATOR TO 

C COMPENSATE FOR THE FACT THAT EQUIVALENCE OOES NOT REORDER COMMON 

C DIMENSION At 3) .0*131 .*131 

CA.l CROSS! A. DA. B I 
P 2 ■ DO T | A , A | 

P > 0SURTIR2I 
V 2 ■ DO T IDA, DAI 
B 2 a OO T t B . B ) 
bB a DSURTIB2) 

A *D0 T I A.DAI/U 
P a B 2/ L 
C 3 « V2-2#0D0*U/R 
SM A *- U/C 3 

ECC aOSURTI 1 .OOOFC 34P/U) 

01 VC a ARKTNSI 18C.BI 31 .D SORT IB II I 442FB{2»*42) 1 
OMGa ARKTNSI 360 .-BI21 ,B II I) 

RCA ap/| 1.OD0FECC I 
THTaARKTNSt 360. IP-R I *BU • A ) 

BET a ARKTNSI 360 • A (2 I *9 I I I -A 1 1 1 *B 12 I , A I 3) *98) 

BE* * BET-THT 
IFIBEP) 2.3.3 

2 BEP*BEPFCa28Jl€9 30 0 

3 CONTINUE 
HTD«5T.i 95779SID0 


1117000 
1 118000 
1119000 
1 120000 
1 121000 
112*000 
I 123000 
I 12*000 
1 12’. 0 90 
1 1 2 DO 10 
I 12 7000 
1 1 2 8000 
1 129000 
1 130000 
1131000 
1132000 
I 133000 
I 1 JAOOO 
1139000 
1 136000 
1 137000 
1 130000 
1 139000 
1 140000 
1 141000 
1 142000 
1142000 
1144000 
I 149000 
1 1 POO 80 


41 


0 INC«OlNCARTO 

OAGMIMGAP TO 

BEP.BEPAPTO 

■»ER V-OSOPTIC J* 2 » OU 0 • 'J/P C A ) 

9t4PaPfcRV-0SURT(U/RCA> 

CTAS-1P/R- 1 .00 01 /ECC 

IMCTAS.LT. UJO.OH. ICTAS-UU0l»Ut.l ( O-6l 60 TO 200 

STAS-O.OC 

63 TO 201 

2C0 STAS«OSQRT( 1 .00 0-C TA SAC TA SI 
STAS«OSICN< STAS. A) 

201 CONTINUE 

THt>ARKTNS< I SO .CTAS.STASI 
THE T» T ME*RTD 

CAwL ICONIC < U.ECC • SMA ,P .THE .TPtR ,F AC» 

TPER>TPf R/6640C.0U0 
IF ( SMA ) If. 10. I I 
10 CONTINUE 
63 TO 12 

It F|bAAO/^SOR1(U*SMA) 

F2« 1*00 0- A /SMA 
SINE “F | /ECC 
C3 St «F 2/ECC 

E>ARKTNS< 350. COSE .SINE! 

X44 9«f E-ECCASlNE IARTD 
12 CONTINUE 

3E( I I • ECC 
3E I 2 1 » SMA 
JF.I 2) ■ OINC 
3E( 41 > OMG 
OEI SI ■ BtP 
3E161 » THET 
RETURN 
END 


141000 

140000 

149000 
150000 
1SIOOO 
ISiO'O 
ISJOOO 
1 34000 
1SS000 

iseooo 

isrooo 

isoooo 

1S9OO0 
1 60006 
161000 
142000 
16JOOO 
164000 
16S000 
166000 
167000 
166000 
166000 
170000 
171000 
I 72000 
1 73000 
174000 
1 76000 
176000 
17/000 
1 76000 
176C00 
160000 


C 


C 


12 


13 

14 


1 1 


16 


SUBROUTINE TCONIC IU .EC .A . SLR .TA 2 .T.FACl 
IM»l ICIT REAL681A-N ,0-21 

TAM 01 QOOCFL » -O SI NIOC OOF L) /DC0S1 QOOOFL) 

AS *046 £1 A I 

FAC -ABAOSum T( 4U/U) 

eca •»( i*rso-tc 1/1 I * 00 OaF C l 

A9E »0 SORT (UAH S( ECA1 ) 

THE «TAN O t • SO 06 TA 2 1 

IF! ABE-.OOCOSOO >11.11.1? 

CONTINUE 

EC A >2. OOC AD A TAN ( ABE A THE I 
IF! A I 14. 1 I . 1 i 
T-FACA( EC A- EC A OSIN(ECA) 1 
CO TO IE 

ANG«. 73 53961 <30 0« • SOOAfeCA 
T *F AC A 1 ECATAN01ECAI — OLOO ( TA NO ( A NO I I I 
00 TO 16 

FAC *OSURT( SLRAA3/1 I A2.000/1 1 1 • 0O0 ACC I • A2 I 
EC 1 *ECA A 1 HE A A 2 

T *PACA ( THE ♦ TMC ** T* t ft. OOO -2 *000 «E C A 1 /3.0 00 - 1 2 • 0 DO -3 .OOO AECA 1 AE .1/5 
l.OOOM 3*0U( 4.0OOAECA) AEC1 A A2/7* 000 -14* 000-5*3 DO AECA) AECIAA3/' .CuO 
A | | 

CONTINUE 

RETURN 

ENO 


101000 
162000 
163000 
164000 
165000 
1660C0 
167000 
160000 
169000 
190000 
191000 
192000 
193000 
194000 
195000 
195000 
197000 
196000 
199000 
200000 
201000 
202000 
203000 
204000 
203000 
206000 
20 7000 


SUBROUTINE HISTO ( T I TLE .NAME ,P • X • AM • XS . AP. XNOM, N3 I 1206000 

OIRCA.SION TITLE (261 .NAME12I .Pll fc.Xl 1 I .AL1 150 I 1209000 

DAT'. BLANK/* •/ 1210000 

OA* A E TE /• I • / 1211000 

•FITE I5.60CI TITLE 1212000 

•RITE 15.6101 NAME. XNOM, XM.XS, NS 1213000 

• RITE (6,620) -SCI I. PHI .Pll 1 1214000 

PKAX * P( 1 I 1215000 

03 10 I * 2.NP 1215000 

PMAX ■ AMAX1 (PMAX.P1I I |T 1217000 

10 C3NT1NUE 12160C0 


DP 

■ *05 








1219000 

IF 

(PMA X*OT *2*5} 

OP 

• 

• 1 





traoooo 

IF 

1 PMA X*OT*5*07 

OP 

8 

• 2 





1221000 

IF 

(PMAX.GT. 10* I 

OP 

m 

• 5 





l 222000 

IF 

<PMAX.GT.rS. r 

OP 

W 

1. 

— 

~ — * — 

— - — — — 


1223000 

IF 

( PMA X.GT.50. 1 

OP 

8 

2. 





1224000 
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03 20 I i 

AIM I ) - 
20 CONTINUE 


> 1*50 

BLANK 


RETURN 

FORMAT ( 1H1.4X.20A4I 


600 FORMAT ( 1HI.4X.20A4I 

620 F 3 R MAT {/5X.» PC Ta • • 50 X • 5 * a9X • • I NTERV ALS • . I 2 » • • PCT • FREQ*'a&«a 
2 * SUM * /5X • * FREU* • .SOXaSXa* I MINUS I hF* • , I PE 1 2* 4, 0P2F I 2 * 7 ) 
610 FORMAT ( 5X. 2A4.2X, *nOMI NAL-* .1 PE12.4 .2X .* ME ANa • , 1PE12.4. 2X . 

9 I < ■ r.Hl a * . IOC I .3 II . I CiUDI Cat .1*1 


2 • SIGMA a • # 1PE I 2.4 .2 X . • SAMPLE - * *161 

650 FORMAT ( 5X.F5*.' *1 X.50AI ,5X,I 3 *1 X *1 P2EI 2*4 ,0P2F12*T ) 

660 FORMAT ( 5X *F 5* 2 • 1 X , 50A1 .5 X* I J « 1 PE 1 3*4 • • PLUS INF* **0P2F12* 

F3RMAT (5X.6X.2015) 

FORMAT (5X.6C. 5X • * 1 N T E R V A S') 

END 


670 

660 


7) 


12X5000 

1220000 

1227000 



A W N a OP* 5C. 



12X6000 


SUM ■ P ( 1 ) 


* ** — • — 

1229006 


03 50 1 • 1 a 50 



12X0000 


03 30 J a 1«NP 



12X1000 


IF (P( J )«DP/2**GE*ALNJ AL1UI a EVE 


1232006 

30 

C3NTINUC 



1233000 


IP 1 • 1*1 



1234000 


SUM a SUM FP (1*1) 




IF ( | *G1 *NP-2) 

GO TO 40 


1236000 


■ RITE (6*6501 

ALN.ALlalPI aX(I),X(|F)| 

. P ( I F I ) • SUM 

• 23 7000 


G3 TC 46 



12X6000 

40 

IF ( | .6T.NP- | > 

GO TO 45 


1239000 


■ titi (e.660> 

ALN .A LI .1 PI .X(l) .P (I F 1 ) 

aSUM 

1240000 


GO H) 48 



1241000 

AN 

■RITE ((a 65C ) 

UN • A L 1 


1242000 

*6 

ALT a ALN-DP 



1X4X000 

50 

C3NTINOE 



1244000 


■ RITE (6.670) 

( I a I * 5 a NP .5 ) 


1245000 


■ RITE (6.680 



1246000 


1247000 
1240000 
1249000 
1250000 
125 1000 
1252000 
1 2! JO 00 
1 254000 
1255000 
1256000 
1257000 


SUBROUTINE OcLVS ( A . R A , RP . A S a t S . A 1 . AI S . APF . AMU . OV S ) 
IMPLICIT REAL*6(A-N .0-2) 

RAS • AS»(I .OOfESI 
RPS > ASX 1.00-ES) 

VAS - O SQR T ( XMU* ( 2* OO/R AS-I . OO / AS ) ) 

VPS > O SQR T ( XMU* (2*0 C/RPS — 1 .00/AS) ) 

ANA a (PA*RPS)/2*D0 
AMP ■ ( RPFRA S) /2.D0 

VA a O SQR T ( XMU* (2* DO/RA— I * DO/A )) 

VP - U SORT! XML* (2*O0/RP-l *00/4 )) 

VAN 1 a O SQR T ( XMU* (2 a 00/RA— 1 a OO/ ANA ) ) 

V*M 1 s DSORT( XML* (2*D0/RP-l *D0/ANP) ) 

VA N 2 > OSQRTI XMUB (2 * DO/RPS- 1 a 00/ ANA I ) 

V* N 2 > 0S0RT( XMU* (2«00/RAS-1*00/ANP) ) 

OVA 1 a OAB S( VAN 1— VA ) 

OVA 2 a DABS! VAN2- VPS) 

OV* 1 a OABSi VPN 1- VP) 

DVP 2 *• UAd S( VPN 2- VAS) 

IF (A I S*GE • 0*00 ) GO TO 100 

OVS ■ OM 1 N 1 (0 VA l*OV42.D VP1 F Q VP2 ). 

RETURN 

IOC CAPF a DCOS(APF) 

R 1 — rp bra /A /( 1 .OOMRA/A -1*00)*CAPF) 

R 2 a RP BRA /A / ( 1 *0 v (RA/A -1*D0)*CAPF) 

R 3 B HP S*R A / ANA / ( i * Du* ( RA /A NA — 1 • OO ) *CAPF ) 

R 4 - RPS*RA/ANA/( 1* OO- ( RA/ANA-1 a DO ) *CAPF) 

R 5 a RP*RAS/ANP/( UOOB (RAS/ANP-1 aOO ) *CAPF) 

R 6 « RP*PAS/ANP/( 1*00- (RAS/ANP-1 *00 ) CCAPF) 

R 7 - RP S*RA£/AS/( laOOB (RAS/AS-1 • 00) *CAPF 1 

R 8 * RP S*RA$/AS/( 1 *00- ( RA S/AS— 1 a DO ) •CAPF ) 

VI a o SORT ( XMU* (2* OO/R 1-1* DO/ A )) 

V 2 a D SQR T ( XMU* ( 2* O 0/R2-1 • DO/A )) 

V 3 a OSQPT ( XMU* (2* 00/R3-1* 00/ANA) ) 

V 4 a OSORT ( XMU* ( 2* 0 0/R4- I a OO/A NA) ) 

V 5 - OSORT ( XMU*(2*D0/R5-1*D0/ANP) ) 

V 6 a OSORT ( XMU* (2*00/R6-l* DO/ANP) ) 

V 7 « OSORT ( XMU* (2.D0/R7-I. 00/AS) ) 

V 8 a OSORT ( XMU* ( 2a Ofl/R 6-1 a 00/AS) ) 

D1 * OAUSI A I-AI S) 

DVS * OM IN 1 <01 «Vl*DVAt*D VA2 * 

2 3I*V2*1>VP|*0VP2 * 

3 0 1 * V5*0 VP1 BO VP2 * 

4 01* V8FO VA 1*0 VA2 a 
RETURN 
END 


D f *V2 fDVAIf0vA 2« 
OI 4V3BOVA1 FOV A2 a 
01 *V6 FOVPl F DVP2 a 
Ol *V7*0VPl *0VP2 a 


Dl*VlFOvPI*OVP2. 
01 *V4 aOV A 1 FOVA2. 
01*V7*CVA1 *0 V a2. 
OI*V8*OVPI*OVP2 ) 


1256000 
1259000 
1259100 
1259200 
1269300 
1259400 
1261000 
1242000 
1263C00 
1264000 
1265000 
1266000 
1267000 
1266000 
1269000 
1270000 
1271000 
1278000 
1273000 
1274000 
1216000 
1276000 
1277000 
1276000 
1279000 
1260000 
1261000 
1262000 
1262100 
1262200 
1263000 
1 28 f ,7 
1266000 
1266000 
1267000 
1266000 
1266100 
1266200 
1269000 
1290000 
1291000 
1292000 
1293000 
1264000 
1295000 
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